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A UNIQUENESS THEOREM FOR TRIGONOMETRIC 
SERIES 


By 8. VERBLUNSKY 
[Received 24 November 1930] 


. THE object of this note is to prove the following theorem: 


THEOREM A. Let the series 


i+ a) 
sag t+ > (a, cosna + b, sin nz), a, =—o0(n),b,=o(n), (1) 
n=1 


have its upper and lower Poisson sums P(x), P(x) finite at every point, 
and let 
P(x) > (x) (2) 
where s(x) is integrable in the sense of Denjoy-Perron in (0, 27). Then 
(1) ts a Fourier-Denjoy series. 
ao 
If P(r, x) = }ay+ > (a, cos nx + 6, sinnx)r”, 
n=1 
then the result remains true when the Poisson sums are not finite at the 
points of an enumerable set E, provided that for x c E, the condition 
lim (l—r) P(r, x) = 0 (3) 
rl 
is satisfied. 

This theorem is related to some investigations given in a previous 
paper,* a knowledge of which will be assumed. The object of that 
paper was to obtain the analogues of certain uniqueness theorems 
due to Zygmund,} when Zygmund’s condition (/'), namely, that 


io 2) . 
> a, cos nx + b,, sin na 
=1 


n 


be the Fourier series of a continuous function, is replaced by the 


conditions a, =0(n), b,, = 0(n). 


In that paper the analogue of one theorem of Zygmund was not 


obtained. This omission is filled by the present note. 
The enunciation of Theorem A renders it necessary to observe that 


* Verblunsky, Proc. London Math. Soc. (2), 31 (1930), 370-86. Referred 


to as I. 
+ Zygmund, Math. Zeit. 25 (1926), 274-90. 
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82 S. VERBLUNSKY 
the method of proof shows that if 4(x) is integrable L, then (1) is 
a Fourier L series. 

We shall suppose in all that follows that a, = 0. 

2. In this section we are concerned with three preliminary 
lemmas. 

Lemma 1. Jf F(x) is integrable L, and 

Fission us S a,, COS wai 3 b,, sin nx (4) 
n* 
is at x = € summable by Poisson’s method, then 
DF(g)>Pé), D®F(€) < P(é). 

This lemma is due to Rajchman and Zygmund.* The enunciation 
is more general than theirs, but this is the result which is established 
by their proof. 

Lemma 2. If F(x) be wpper semi-continuous in an interval (x, B) 
containing the interval (a, b) in its interior, and if D?F (x) >0 at every 
interior point of (a, B), then 

r—a, 


F(x) < F(a)+*—“(F()—F(@}  (a<x <b). (: 


' b—a 


CU 
~— 


If not, there is a point ¢ such that 
F(c) > F(a)-+—“{F(b)— F(a)} (a<c<b). 
h—a 

Writing O(x) = F(x)+e(x—a)? (e > 0), 


we can choose « such that 


O(c) > O(a) +°—“{@(b)—(a)}. 


gf 


¢— 


Consider the function 
b(x) = (x) —O(a)— ((b) —O(a)}. 
We have h(a) = $(b) = 0, d(c) > 0. 
Further, ¢(2) is upper semi-continuous. Hence it attains its maxi- 
mum at a point é satisfying a<& <b. At &, 
D*O(£) = D*4(é) < 0. 

But D°@(é) = D?F(é)+-2¢. 
Hence D?F(é) < —2e, 
which contradicts the hypothesis. 

* Rajchman and Zygmund, Math. Zeit. 25 (1926), 261-73. Cf. in particular 
the footnote to their Lemma II. 
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Lemma 3. If F(x) be upper semi-continuous in an interval (a, B), 
and if D? F(x) > 0 at each point of the interval, then F(x) is continuous 
in the open interval (cx, B). 
From the relation (5) we infer that, for «a<#<y<z<f, 
F(x)—Fly) — F(x)—F() 
= 
x—y xr—z 
or, what is the same, 
F(x)—Fly) — F(y)—F() 
enn ' 
x—y y—z 
Hence, fora<x<y<z<w<f, we have 
F(x)—Fly) — F(z)—F(w) 
< : 
x—y z—w 
Let 2 tend to y and z to w in a suitable manner and we infer that 
DF (y) < DF w). 
In other words, D F(x) is a monotone-increasing function in the open 
interval (a,f8). A similar argument shows that each of the other 
extreme derivatives is monotone-increasing. Hence, in any interval 
strictly interior to («,8), all four derivatives are bounded. This 
proves the lemma. 
3. We turn now to the proof of Theorem A, and consider the case 
in which the Poisson sums are everywhere finite. We may write 


2) 


in — b | a, COS NX + b,, ain ne (6) 
= n? 

the series being everywhere convergent. Then F(x) possesses the 
property R (I, Lemma 4): i.e. the points of discontinuity of F(x), 
with respect to any perfect set P, are non-dense on P. Further, by 
Lemma l, D2 F(x) > P(x) > h(x) (7) 
at every point. 

We consider an arbitrarily assigned interval A = (a,b), and sup- 
posing, as we may, that ¢(x) is everywhere finite, we denote by ¢() 
a minor function of (x) in (a,b). Writing 


(x) = J o(t) dt 


and G(x) = F(x)—®(z), 


we have D?G(x) >0 (8) 
G2 

















84 S. VERBLUNSKY 
at every interior point of A. Further, G(x) possesses the property 
R within (a,b), and is, like F(x), a differential coefficient. 

We now show that G(x) cannot have an isolated discontinuity 
within A. For suppose that é is such a discontinuity. Then there is 
an «€ >0 such that € is the only point of discontinuity of G(x) in 
(Eé—e, +e). Since G(x) is continuous in ({—e < x < &), it follows by 
the proof of Lemma 3 that G(x) has increasing derivatives in that 
interval, and hence, that there is a constant K such that 


G(x)+ Kau 
is monotone-increasing in (€—}e < x < &). 
In a similar manner, we can find a constant ZL such that 
' G(a)+ La i 
is monotone-diminishing in (€ << # < &+-4e). 
The function H(x) defined by 


H(a Har) + (€—ke <a <6) 
H(é) = G(é)+-Ké 
H (ax) = G(x)+La+(K—L)é (E<a<é+he) 


is continuous in (€—}e <a<é+}e). From the definition of this 
function it follows that, G(a) can have only an ordinary discontinuity 
at €, and since G(x) is a differential coefficient this is impossible. 

4. We have thus proved that the points of discontinuity of G(x) 
on A form a set E which is dense in itself; and this set # is every- 
where dense on a perfect set P. 

We can in virtue of the property F find an interval dc A the end 
points of which are points of P, such that G(x) is continuous on Pd, 
while its points of discontinuity on d are everywhere dense in Pd. 

Let («,,,8,,) cd be an interval contiguous to Pd. In the open interval 
(x, B,) G(x) is continuous, and by (8) and the proof of Lemma 3, 
has increasing derivatives. It follows that G(x) cannot have a dis- 
continuity of the second kind at a, on the right, nor at £8, on the 
left. As a differential coefficient it cannot have an ordinary discon- 
tinuity, so that G(x) is continuous in the closed interval (q,, 8,,). 
Furthermore, from the fact that G(a) has increasing derivatives 
within («,, 8,,), we infer that G(x) must attain its upper bound in 
the closed interval («,, 8,,) at an end point. Since, then, G(x) is con- 
tinuous on Pd, it must be upper semi-continuous on d. By (8) and 
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Lemma 3, it must be continuous on d. This contradicts the assump- 
tion that # exists. Hence G is continuous in A. 

In virtue of (8) and the proof of Lemma 3, 


G(x) = F(x)— | d(t) dt 


has increasing derivatives within A. If we choose a sequence of minor 
functions ¢,(x) which tend uniformly to 


[ ye dt, 


ee 
we infer that F(x)— | dy ( y(t) dt (9) 
a a 
has increasing derivatives within A. It follows that the expression 
(9) has almost everywhere in A a second differential coefficient which 
is integrable L in (a+e, b—e). But 
oa y 
[ dy { ye) ae (10) 
a a 


has almost everywhere in A a second differential coefficient (x) 
which is integrable D in A. Hence F(x) has almost everywhere in 
A a second differential coefficient which is integrable D in (a+, b—e). 
Since A may be arbitrarily chosen, it follows that D?F(x) exists 
almost everywhere and is integrable D in every finite interval. At 
a point € at which D? F(x) does not exist, there is, by Lemma 1, since 
P(é), P(é) are both finite, a finite number between D?F(é) and 
D®F(é). Hence, by a classical theorem of de la Vallée Poussin,* (1) is 
a Fourier-Denjoy series. 

5. We consider now the case in which the Poisson sums are not 
both finite at the points of an enumerable set HZ, at which the con- 
dition (3) is satisfied. This requires a somewhat difficult argument. 
We begin by stating the following lemma which is implied by the 
work of Zygmund (loc. cit., Theorem ITI). 

Lemna 4. Under the conditions of Theorem A, if (x, B) be interior to 


an interval in which 


a, cosnx + b, sin nx 
re Zz n> 
converges to a continuous function F(a), then D® F(x) exists almost 
everywhere in (x, B) and is integrable D in this interval. 


* de la Vallée Poussin, Bull. de ! Acad. Roy. de Belg. (1912), pp. 701-7. 
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We next prove 
Lemna 5. If at the points of an interval (a,b) containing the inter- 
val (a, B) in its interior, the upper and lower Poisson sums of (1) are 
finite except at the points of an enumerable set E; if, for x not belonging 
to E, 





where us(a) is integrable D in (a, b), while for x c E condition (3) is satis- 
fied; if further the series 


F(x) = — > a, cos nx + b, sin nx 


2 


n* 


converges everywhere in (a, b), then F(x) is continuous in (a, B). And, 
if d(x) denote a minor function of (x) in (a, b) and 


(x) = | d(t) dt, 


on D*{F(x)—(x)} > 0 
at all points of (x, B). 

The set E of points at which the Poisson sums are not both finite 
is an ordinary inner limiting set.* Being enumerable, it contains 
no component which is dense in itself. We can write 


E=F'= £,+£E 


where ZL, is the set of isolated points of H! and is not null. Similarly, 
we can write, if EH? is not null, 


BE? = BE, BS 


where EF, is the set of isolated points of H? and is not null. Pro- 


ceeding in this manner, we have since EZ is enumerable, 
E= > Ez 

— 

B<y 
where y is some ordinal of the first or second class. 

Consider a point x,c #,. There is an e > 0 such that (%)—2e, x)+ 2c) 

contains no point of # other than x». For all x in 6 = (x—e, %)+e) 
we have, by Lemma 1, 


D? F(x) > (x). 
Hence. D»* F(x)—®(x)} > 0 (2 #2, 2% <4). (11) 


* Cf. Zygmund, loc. cit., footnote, p. 285. 
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As in § 3, 4, we can show that 


G(x) = F(x)—(2) 


is continuous in (%—e <x <4) and in (%<#<a,+e). Further, 
by the argument at the end of § 3, G(x) is either continuous at x, or 
has an ordinary discontinuity at x). The latter is impossible. Thus 
F(x) is continuous in 8. 

We now apply Lemma 4 to show that (11) is satisfied at 2. By 
that lemma, D?F(x)= g(x) exists almost everywhere in 8 and is 
integrable D. In (x»—e <2 <q) we have 


z= y 

F(x) = f dy i g(t) dt+Ax+ B, (12) 
and in (% <x <2 +e) we have 

F(x)= | dy [ g(t) dt+-A’x+ B’. 
x y 
—_= — 
Hence lim = h)—2F (x9) 
h—0 


exists. It can easily be shown that the existence of /(x,) implies that 


lim (l—r) P(r, 2%) = U(29). 


r—>1 





= (x9) 


By (3) we have /(x)) = 0. It now follows that (12) holds for all x in 
. F(x-+-h)+F(e—h)—2F (x) 


lim 


h—0 h 


= 0 





for all x in 6. Clearly 


lim Pe L) TO —A)— 20 (2) _ 9 
h—0 h 





for all x in 8, so that G(x) has the same property. If, then, (11) did 
not hold at x», it would not hold at an unenumerable set of points 
in 6, a conclusion we know to be false. 

It thus follows that F(x) is continuous in any interval which con- 
tains no point of H?, and that (11) holds at every point of that 
interval. We consider next an isolated point of H?, say xc H,. The 
above argument can be repeated precisely. By transfinite induction 
we thus establish the lemma. 

6. We now show that under the conditions of Theorem A, the 
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series obtained by two formal integrations of (1) converges every- 
where. We first observe that the continuous function 

a a, sinnx — b, cosnx 

B(e)=— >, = (13) 


has a differential coefficient equal to 
a S a, cos nx + 6, sin nx (14) 
n2 
for x not belonging to HZ. We can show that for éc EH, the four 
derivatives of §(x) are bounded. For in virtue of (3) we have 


3 a, cos ng + b, sin ng inca 0 (log 1 ) 
n . Ir, 


and 


2 


FY 2008 né + 6, sin n€ on ~ 0(1). 
hos n 
Since the coefficients in the last series are o(1/n), the partial sums of 

a,, cos n€ + b, sin n€ 

3 7 n* 

are bounded. It follows* that the limits of 

B(E+h)—H(E—h) 

2h 

as h->0 are bounded. On the other hand 


_ B(E+h)+B(E—h)—28(€) _ 0 
h—0 h 
since the coefficients in (13) are o(1/n”). This establishes our assertion. 
We denote the sum of (14) at the points at which it converges by 
F(x). We proceed to show that it converges everywhere. The set of 


points of non-convergence are contained in E and therefore contain 
no component which is dense in itself. This set therefore has an 
isolated point €. We can find an « > 0 such that € is the only point 
of non-convergence in (€—2e, €+2«). Let 6 =(é—e, +e), and let 


¢(x) be a minor function of ¢(x) in 6. Then by Lemma 5, F(x) is 
continuous in (€—e < a < €) and in (€ <<a < +e), and 
D*{ F(x) —®(a)} > 0 (a A€, ac 8). 
Writing G(x) = F(x)—Q(z), 
it follows, as in § 3, that there is a constant K such that 
G(x)+ Ka 

* Cf. the argument on p. 636 of Hobson, Theory of Functions of a Real 

Variable, vol. ii. 
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is monotone-increasing in ((—« <x < €). Hence G(é—0) has a mean- 
ing, and is finite or infinite. The same applies to F(€—0) and to 
F(€+-0). Since the derivatives of §(a) at € are bounded, F(é—0) and 
F(€+-0) are both finite. Hence (14) is summable by Poisson’s method 
at x = &, and so, by Tauber’s theorem, is convergent at that point. 

Since (14) converges everywhere, denoting its sum by F(z), it 
follows by Lemma 5 that F(x) is continuous everywhere, and so by 
Lemma 4 that D*F'(x) exists almost everywhere and is integrable D 
in every finite interval. Writing D?F (x)= f(x) wherever the first 
member exists, we have, in any interval free from points of EZ, 

x y 
F(x) = [ dy | fc) dt+-Ax+B (a<a<b). 


a a 


By the argument on pp. 382-3 of I, we can show that this relation 
holds in every finite interval, and this establishes the theorem. 
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A PROBLEM CONCERNING THE HYPER- 
GEOMETRIC EQUATION (II*) 
By E. G. C. POOLE 

[Received 13 January 1931] 
7. ResumineG the investigation of cases where the hypergeometric 
equation has n distinct solutions, whose product is monodromic, we 
consider in this paper the values n = 4 and n= 5. It appears that, 
for all values of n, there is a solution of the problem in which the 
equation has an apparent singularity and admits two independent 
monodromic integrals, and another solution when the equation has 
two quadratic singularities and is reducible to the form considered 
by Elliot and Darboux. Apart from these solutions, which are 
explained in (8) and (9), and which were already prominent in the 

previous paper, there is no solution for n = 2, 3, or 5. 

For n= 4, there is a different kind of solution, when the para- 
meters of the hypergeometric equation differ by integers from those 
of certain reduced forms. Two singularities at least must be of the 
third order, the remaining one being of either the second or the third 
order. But this condition, though necessary, is not sufficient. By 
elementary processes, taking the points « = 0 and a = 1 as points of 
the third order, and assigning the exponent zero to the appropriate 
solution at each of these points, the reduced forms, in which we take 
l>a>b>0, 1>c>0, are found to be 

F(3,3;35%), = F(3,3595%), = F(4,35$52), FF (G.45 952). 
The examination of this solution occupies the greater part of (10); 
the remaining paragraphs of (10) and the whole of (11) are occupied 
in reviewing and rejecting all other possible types of solution. 

8. It was seen in (3.2) and (5.4) that, when one singularity is 
apparent, the other two are of the same order, and the type of the 
equation is (V,N,1); excluding equations with logarithmic points, 
which cannot yield a solution of our problem, there are two indepen- 
dent monodromic solutions, which may be taken in the canonical form 

y, = F(a,b;c;2), Yo = x1-*(1—2x)°-* F(1—a, 1—b; 2—e; 2), 
* Quart. J. of Math. (Oxford), 1 (1930), 108-15. We take the opportunity 


of noting the following corrections: p. 111, line 4, for (2.6) read (2.16); p. 114, 
line 3, for g?= 1 read s§=1; p. 114, line 5, for g=1 read s=1. 
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where a, b are integers, a << 0 << 1<), and ¢ is rational, its denomina- 
tor being N, the order of the singularities. Conversely, two mono- 
dromic solutions imply an apparent singularity. 

It is clear that (Ay’+ By’) is a monodromic product, and is an 
irreducible solution of our problem for n = N. Expressions such as 
y,(Ay + By) or y, y.(Ay* + By’) are also monodromic solutions 
for n = N+1 and n= N+ 2 respectively. But as they break up into 
distinct groups of factors, which are separately monodromic, they 
may be termed reducible, and will be left aside as trivial. 

9.1. When there are two quadratic singularities, it was seen in 
(4.2) that, instead of Gauss’s form, we could use the equation of 
Elliot and Darboux, 


- = [m(m-+-1) cosec?@-+-n(n-+- 1) sec?6+-k*]y, 
where the independent variable of Gauss’s equation is 2 = sin*@, and 
the exponent differences are m+4, n+}, and ik, m and n being 
‘integers. Thus 7k is a real fraction, whose denominator N is the 
order of the third singularity in the symbol (2, 2, V). 

The solutions can be written 


y, = & R(cot 6), Yo = e-** R(—cot 8), 


where the symbol FR represents a rational function, the only poles of 
R(z) being z = 0 and z = oo. 

The conditions for monodromy at the points x= 0, x= oo, and 
« = 1 of Gauss’s form, correspond to the conditions that our product 
(which is a homogeneous polynomial in y,, ¥2) shall be monodromic, 
when @ is replaced by —@ or by 6+-7. 

The polynomial must accordingly be monodromic, when (y,, ¥.) are 
replaced by (Yo,y,) or by (wy,,wy2), where w = e*7 and w* = —1. 
After verifying that the polynomial satisfies the conditions, we 
examine whether it breaks up into factors, which are separately 
monodromic. Such solutions would be reducible and we have agreed 
to ignore them. We tabulate the irreducible solutions, leaving the 
verification to the reader. 

9.2. If n=4, we have N<4 by (5.3). The first substitution 
separates the two types of solution 


A(yj+y3)+ By, yolyi+y3)+ Cyj yi, 
(y7—y3)[A(y7 +93) + By, yo). 
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The second substitution yields, after eliminating the reducible cases: 
N=1, A(y{+y$)+ By, yo(yi+y3)+ Cyj y3, 
WY = 2, 3, no irreducible solution, 
N=4, A(yj+y3), A(yi—y?)- 
9.3. If n= 5, we must have N < 5. 
The first substitution separates the solutions 
A(yi+-y3)+ By, yolyi +y2)+ CyZ Y3(Y1t Yo); 
A(yi—y2) + By, yolyi—Y2) + CY{ Y3(Y1—Y2)- 
The second substitution yields no irreducible solutions for N = 1, 2 
3, 4. If N = 5, the irreducible solutions are 
A(yj+y3), A(yi—y%)- 
9.4. The discussion of this and the preceding section has disposed 
of the cases, where the symbol (N,, N,, N;) contains unity, or where 
2 is repeated. We need only attend to those symbols with no figure 


less than 2, and with two figures not less than 3. Any admissible 


equation falls under one or more of the following types, which will 


now be dealt with: 

n=4,  (4,4,N), (4,3,N),  (3,3,N); 

n=5, (5,5,N), (5,4,N),  (5,3,N), 
(4,4,N),  (4,3,N),  (3,3,N). 

10.1. It was shown in (5.2) that, if y,, y. are the fundamental solu- 
tions at a singularity of order N, the factors y,, y, may occur singly 
in our monodromic product, but all other factors occur in groups 
(Ay) + By). Absorbing the coefficients, real or complex, into the 
solutions, the canonical forms for n = 4 are (Pj, +P},) if N = 4, and 
P,, (Pi, +P2) if N=3. Thus a solution of the type (4,4,N) or 
(4,3,.V) would lead to one or other of the identities 

+L = ite 1h, 2... +h 


: Xe 1 


~ 


respectively. If we proceed as in (6.2), putting 
Py, = pP,, +4P,,. P,,=1P,,+8P,,, ps—qr ~ 0, 
and comparing coefficients, it is found that the second identity leads 
to no solution, and that the first can be reduced, by proper choice of 
notation, to the case P,, =P, P,,=P,,, which is covered by the 
discussion in (8). 

10.2. In (3,3,.V), we use Gauss’s form and we observe that a, b, c 
are real, because the exponent-differences are real. There is no loss 
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of generality in assuming that the singularities of the third order are 
at x= 0 and x= 1, and we can arrange that the branches | | 
Fas Fa shall be real, if x is real and 0<2<1. We have, in that 
interval, the identity 

P, (AP? + BPE) = P,,(CP),+ DP), 
whose coefficients might perhaps be complex. Now, by (2.3), P,, is 
identifiable with one of the linear factors on the right, and for an 
irreducible solution, this must not be P,. Hence there must be 
another real linear factor, so that C/D is real, and similarly A/B is 
real. Absorbing a real coefficient into each fundamental solution, 


we find 3_| D3\—@ 3 1 ps 

P,,(P2+P3) =OP,,(P3, + P3). 
The coefficient © is the ratio of two real functions, and is therefore 
real, and may be taken as +1, without loss of generality. We now 


have P,,=UP,,+P,), P, =m(P,,+P,,), (10.21) 
where /, m are real and different from zero. Hence 
he —F Fie) = Lee FF 
B(3m?2P2 —3mP, P, +P?) = +m*(PE,—3lP,,P,,+3?P2 ). 

The coefficients of P,,P,, give ? = +m?*, but the lower sign leads to 
imaginary values, and must be rejected. Our identity now requires 
2? = m?, l= +m, 3lm = 1, 
but again the lower sign leads to imaginary values, so that we have 
l=m= +1/N3. 

These two forms differ only in notation, so we may take the upper 

sign. We must now verify that the formulae 

V8P,=P,+P,, 8P,,=P,,+Py (10.22) 
are compatible with the classical theory of the hypergeometric equa- 
tion. On removing the factor x*(1—x)”1, the branches P,, and P,, 
are identified with those of exponent zero at «=0 and at x= 1 
respectively. We consider the real interval 0 <x <1, and, to keep 
the formulae as simple as possible, we use the notation :* 


P _xY, = F(a,b;c; 2), 


a 


P,, ©Y, = 2!~F(a—c+1,b—c+1;2—c;2), 


P,, <Y;= F(a,b;a+b—c+1;1—2), 
P,, © ¥, = (1—2z)** °F (c—a, c—b;c—a—b+1; 1—2). 


(10.23) 


* Cf. Forsyth, Differential Equations, Ch. VI. 
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Now the formulae (10.22) must be equivalent to 


I'(c) \(e—a—b) (c)'(a+b—c) 


, ( a 

7 l(c—a)T(e =p ot 1) Vy 
(1 
(1 


(a) (6) 
M(i—e)P(+a+b— ¢) y I'(c—1)l'(a+b—c+1) 
l(1+a—c)l(1+b—c) T(a)(b) . 
On comparing these two sets of formulae, we find: 


V3P,, _T'(c—a) T'(c—6) Y, V3P,  I'(1+a—c)I'(1+b—c) Y, 


i= ¥,+ 


y. a6 Tc—a—b) ¥,’ P,, Ti—e)PU+a+b—e) Y,’ 
hence 
_ Te—a)P(1+a—e)P(c—b)P(1+6—c) _ sinzesina(c—a—b) 
(ce) (1—ce)T(e —a—b)T(1-+-a+b—c) ~~ sin 7(e—a)sin 7(e—b) ° 
(10.25) 


In this equation, we observe that, if a, b, c is any solution, another 
can be found by adding any integers to a, b, c. We shall therefore 
fix our attention on the solutions for which a, b, c are positive proper 
fractions, which will be called the reduced solutions. Since a, b occur 
symmetrically, we may take 
l>aS>é >, l>c3>0. 

We may note that the hypergeometric function with a=b would 
itself be logarithmic, and so could not lead to a solution of our 
problem; but an equation whose reduced form was of this type might 
have an apparent singularity. Actually, the case does not arise, as 
(a—b) is found not to be an integer. Now, since x = 0 and x = 1 are 
singularities of the third order, we must have 


See » 2 — 1 
c= or §, c—a—b = p+}, 


where p is an integer. If these expressions have the same fractional 
part, so that (a+b) is an integer, (10.25) gives 


3 sin 7(c—a)sin 7(c—b) = cos 7(a+b)sin?xc = }cosm(a+b), 
2 cos 7(a—b)—2 cos n(2c—a—b) = cosx(a+b), 
2 cos 7(a—b) = cos 7(a+-b)| 1+ 2cos 27c| = 0. 
The solution requires, with our conventions, 
a—b=}, a+b=1, c=tor &, 
and the corresponding reduced forms eis to 


F(3,4;452) or F(3, 45452). 
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If (a+-6) is not an integer, (2c-—a—b) is an integer, and we can then 
put (10.25) in the form 

3 sin 7(c—a)sin 7(c—b) = —sin*rc cos 7(2c—a—b) 
= — cos 7(2c—a—b), 
2 cos 7(a—b) = cos m(2c—a—b) = +1. 
Taking the upper sign, our conventions give the reduced solution 
a—b=}, a+b = 2c, c=} or , 
and the hypergeometric functions are derivable from 
F(, 33452) or F(%, 3; §;2). (10.27) 
If we choose the lower sign, (22-—a—b) = +1, and (a—b) = 3, and 
there are no reduced solutions. 

The exponent-difference at x= 00 being (a—b), we see that this 
point is of the second order in (10.26), and of the third order in 
(10.27); the symbols of the families are (3,3,2) and (3,3,3) re- 
spectively, although not every hypergeometric equation with these 
symbols belongs to our families. For example, if a = 4, b= 0, c = 3, 
all the exponent-differences are 4, but the equation is not admissible. 

If a, b, c differ by integers from a set of reduced values, we shall 
have, in the region where |x| <1, |l—az| <1, 

oe 3 £5 on 31 ps 
F(x) = PL (Pot Fa) =f, Fy, + *),)- 

‘At x=0 and «= 1, F(x) remains uniform, and is either regular or 
has poles at those points, according as the integers 3(1—c) and 
3(a+-b—c) are positive or negative. It has no other singularities in 
the finite part of the plane, and so must be a uniform analytic func- 
tion which has at most a pole at infinity; it is therefore a rational 
function of x. Equivalent forms of F(a) include 

VSP Fe Fak atta JOP FY (Pe —P Fea 

3V3P, P, (PG,—V3P,,P,, +P), 3V8P,, P, (Fo, — Py) 

Sv Fe. lat we) 

[f we pass from these to the classical hypergeometric functions, we 
can find a number of equivalent forms, which are multiples of F(x), 
such as 
ry (\(1—c)'(1+a+b—c) ,, ~~, I 
Y UREA S ae 
FY AP(1+e—e)I(14+6—e) * am 


which is a rational function of 2. 


“(c)\(c—a—b) v3) 


(c—a)I'(ec—b) * 
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11.1. We now show that none of the six types of hypergeometric 
equation enumerated in (9.4) can yield a solution of our problem for 
n= 5. For the first three types, a formal argument following the 
lines of (6.2) or (10.1) will show that the identities 

(5, 5, N) 1+ =f.+F. 
7 5} _ 4 H 
(5,4,N) P5+P3=P,(P!+P), 
~ 9 W 5 ae 3) D3 
(5, 3, V) the. = P,P, Py, +P), 
yield only the type (5, 5, 1), covered by (8). 

11.2. In the other cases, the reader will argue, as in (10.2), that, 
in the real interval 0 <a <1, we have real identities connecting the 
four real solutions P,,, P,,, P,,, P,,. The proof then proceeds as 
follows, in the three cases. 

11.3. For (4,4, V), the real identity is 

n ie { n 
FLAG Fe) BP ee 
For an irreducible solution, the three real factors must correspond 
in a way, which (except in notation) can only be as follows: 


Py = UP, AP), Py, = mPa, t+ Pa); (ee OS a 


Y a 


We find 4/m = 1, so that we put p= 2/ = 1/2m, and obtain 
2P,, = p(P,, +P,,)» 2P,,, = p(3P,,—F,,); 


p°(P? +P2)-+(P2 +P?) =0, 


Y2 
p°(5P? —2P, P,,+P7,) + 2(P7, +P; ) = 0, 
which is impossible, when P,,, P,, are independent. 
11.4. For (4,3, .V) the real identity is 
4 4 < ° 
Py, (Po, — Pag) P, Py (Py,+Py,): 
and the correspondence of factors, for an irreducible solution, can 
always be written 
e = U(P,., ! Pd = m(P.,. —Fa)> Pix (PY, +P). 
We have / = m, and reduce the identity to 
(P32 +-P2.) = 2PPF,—P,,P, + PF.) = Prat 3P? ), 


which is impossible, when P,,, P,, are independent. 


11.5. For (3,3, N), the real identity is 
Py Pag( Po, +Pe,) = P,P, (Py, +Py,)s 


xy" Og 


and an irreducible solution would require that and P,, should 


which is 


Oy 
be proportional to two distinct real factors of P},+-P},, 
impossible. 





A PROBLEM IN ADDITIVE ARITHMETIC 
By LEONARD CARLITZ 
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Introduction 
Tus paper is concerned with the following problem. 
Let a(n) be a function of the positive integral variable n satisfying 
the following conditions :* 
afl) == I, a(n) = O(né), a(p) = 14+ O(p-**), (1) 
a(mn) = a(m)a(n) for (m,n) = 1, 
where p is any prime, and « is an arbitrary positive quantity. Then 
we seek an asymptotic expression for 
R,(n) = ¥ a(n,)a(Ng)... a(n,) (v > 3), 
the sum being taken over all the partitions of n into v positive integers. 
The method used is that of Hardy and Littlewood, particularly 
that as applied to the Goldbach problem, and simplified by Landau 
in his Vorlesungen iiber Zahlentheorie.t It should’ be remarked, 
however, that it is here not necessary to introduce any assumptions 
relative to the zeros of the Dirichlet L-functions. Our main result 


is the following: 


| Bom) — dil S(n)| < An’-*e, 


| (1)! 
where b= LT (s Bee ie ts a(p* )—o(P) o) 


Pp 
S(n) = an oP) x 
\ Pp 


—p! 
“oe p)+ GP (pl) 


pity 


p p'—p'- 
1+ 
CD T(p)+F er pi 
, (3) 
7’(P) 
p” 
p! being the highest power of p dividing n; 





pin [=_ 


* Slightly less stringent conditions on a(n) would suffice. Thus the condi- 
tion on p need not hold for a finite number of primes; also the 1/2 in the 
exponent might be decreased somewhat. 

+ ‘Some Problems of Partitio Numerorum, III: On the Expression of a 
Number as a Sum of Primes’: Acta Math. 44 (1923), 1-70. 

t Bd. i, 5. Teil, pp. 183-234. Cited as Vorlesungen, i. 

3695.2 H 
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o(p°+!) —a(p? 

a(p*)—atpe-2) + Wa") 
P 


7(p*) = 1, 


Pp 

and A is a positive quantity depending only on e. 

If we suppose that, in addition to the restrictions in (1), a(n) is 
always either zero or one, then (2) gives an asymptotic expression 
for the number of representations of a large integer n by v numbers 
n; for each of which a(n;) = 1. Such an application will be considered 
in § 5. We also give one other application of the general formula (2) 
(§ 6); it is easy to see how many other applications may be made. 








1. Transformation of the generating function 
We are going to follow rather closely the analysis in Part V of 
Vorlesungen, i. We begin with the function 


f2)= S (1m Jac”, (5) 


which is evidently convergent inside the unit circle of the a-plane. 
If G denote a circle with centre at the origin and of radius e~”, then, 
of course, 1 f(x) 
R,(n) = — | —; aa. (6) 
2mt J a 
“ 
Consider the Farey* series //k, 
0<k<.n'i, 0<1<k, @,&) = 3. 
The circle G is divided into ares B, by means of the points 
exp(—*+ rwaand 
ae k+k’ 
1/k and l'/k’ being two adjacent fractions in the series. Each arc con- 
tains just one point e~!"p (p = e”*"/*): on B, put 


, = ] : 
z=pd, a = ¢-%, y= -——h, 


n 
so that * nr 
—0,<0<4,, oe te Oy KE ae i= 2. 2). 
NN a f) 2) 
We shall replace f(x) on B, by 
b 
f(x) = iy") 


* Cf. Vorlesungen, i, p. 217. 
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where 7(k) is defined by (4) and 
t(k) = [] r(p*), io II pr. 
Pe pie 
Our object is to prove 


f(x) —,(x)| < An?8+€. 
By Mellin’s Formula,* for poten >0 


ev = .. : [ y*I'(s) ds, 


ar? 
(p) 
where _ Bt ei 
j= | (B real 
(p) B—wi 
Now P mo—ym 
f(x) = > a(m)p™e- 
1 
=s5—5 (m adit (my)-*1(s) ds. 
2. 
1 a(m)p™ 
= fro SS) 
271 ms 
(2 
But os 
> a(m)p™ => +4 (m)p™ 
ms 
m=1 a\k (m,k)=d 
l x(dm)p dm 
k al (m,d)=1 
1  & x(dm)p (8) 
Kis? (m,8)=1 al 
o(5) denoting a primitive 5th root of unity. If we put 
5. 
B(x) = 2 x(m)p™, 
m= 


y being a character (mod 6), the sum in (9) becomes 


da d* wD B(x >) si ad 


k=dd m=l1 


* Mellin, ‘Abriss einer einheitlichen Theorie der Gamma- 


geometrischen Functionen’: Math. Annalen, 68 (1910), 305-37. 
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(9) 


(10) 


und hyper- 


+ The interchange of integration and summation may be justified without 
g } 


difficulty ; ef. Vorlesungen, i, p. 213. 
ifficulty ; ¢ orlesungen, i, p i” 610915 A 
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Now letting d = Iq’, 


a(dm)x(m) 
2 “ Sy 
-T1 (1+ x( p P)x(P) P) , APP p*) ro) TT ( (q) +e @ ,...) 


78 
pid * I q\d 


=TI(+ sex?) +.) x 


allp *‘ 


qid * 
= n(s, x) €(8, xX: d), say. (11) 


As is customary, put 


becuadl ; x(m) __ lI ( + oy | 


m=1 m p Pr 
then : 
n(8,x) __ - , a(p)x(p) , ) 
#53 ~ T1(+seRe+-)(-92) 
| | (1+ ei x) = = P) (2 >) ) 


n(8, x) a (1 ey ie i? 
L(s, x) 7 I] T Ben - p?9(1—p-*), 
<A = A(e) for o> $+e. (12) 


We see, then, that in the half-plane o > 4, (s, x) is regular except 
for a pole of order one at s = 1 when x = xo, the principal character. 
The residue at this pole is by (12) 


— \¥ a(p)—l1 , bd(d ; 
Ce caren ee 
where a a(p)—2 
we) =| [(1+ ; +~). 


pid 


Returning to (8), we deduce from 


a(m)p™ __ -1) 
ee dre 5 2, Bex me, xi(exd) (24) 
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that the integrand has a pole of order one at s = 1; using (13), the 
residue at this _ is 


 59(9) i. 
= 5 ab aie) PO say fll x04) (15) 


1 £(1, x9, d) * 
be: @ 


ETE ro? 


pik pl=dd 
lec 7) €(1, xo(P), P”™)) 
=m LT {6 E( » Xo(1 ),p’)— b(p) | 


( i 
Poe : —a(p!-1) 


bei (p*) 1 MP)—1 
3 Ba pai 
b 


= Fay (4) = Hole) (16) 





-2U 


Xeturning again to (8), we consider the integral of 


m 
rere) > me 
m* 


taken in the positive direction around the rectangle of vertices 
2+ Mi, d+e+ Mi, (M > 0). 
From the results above the integrand is uniform within and on the 
rectangle and is regular except for the pole at s = 1 with residue as 
given by (16). Furthermore, exactly as in Vorlesungen, i, p. 215, it 
may be shown that the integral along o = }+- converges and that 
the integrals along the horizontal lines approach zero when M be- 
comes infinite. Accordingly 
L[=4+s; | (17) 
2ri J =? 2Qrt J 
(2) (4+) 
We must now find an upper bound for the integral in the right 


member of (17); this requires a lemma on the magnitude of L(s, x) 
on the line o= $+. 


* (8) is the Mobius p»-function ; it appears here as B(x). 
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2. Magnitude of L(s, x) in the critical strip 
Lemma. On o = }-+.¢, (e >0; x, modhk), 
L(s, x)| < Ak*+€(|t|+-1)#+€, (18) 

Proof. We may suppose x a primitive character (mod k). For if 
it is imprimitive, there exists* a K, K\k, and a primitive character 
X (mod K) such that 

L(s, x) = [10- =(P)) L(s, X}. 
\ P| 
Then on o = 3+, 
L(s, x)| s one a 


\L(s, X)| T] 2 < Ak«|L(s, X)}. 
pik 


Hence, by properly modifying ¢«, if (18) be proved for all primitive 
characters it will hold generally. 
If now, x being primitive, we let 


oa s+a 


Bs, x)= (7) @ r(* 7) Hex) 


where a = 0 for x(—1) = 1, a= 1 for x(—1) = —1, thent 


&(s, x) = H(x)=(1—s, x-), H(x)| = 1. (19) 
Let s = —e-+4ti, so that 
| L(1—s, x) | < A(e). (20) 
From (19), 
w\-St¢ (sta ie ‘a\ —! ste _(l—s+a — 
(7) 2 r , )L@,x)=H00(Z) 2 r : )LO—s, x ). 


But 


— ; > eee “\r(A) 





r(° as | r(-3 “yr *) 


1+€+a_} 2+ 


<e™"A(eje Pi |t/ 41-1) 2 *A(e)e-2"(|t|-+1) 
< A(|t}+1)?+¢. 


Hence L(8,x)< 


as . * a(jt|41)* 


\7, 
< Aki+«(\t|+-1)!+€, 
* Cf. Landau, Handbuch der Lehre von der Verteilung der Primzahlen, i 
(1909), pp. 478 ff. 
PI 
{+ Landau, loc. cit., p. 497. t Cf. Vorlesungen, i, p. 195. 
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We now apply an extension of the Phragmén-Lindeléf Theorem, 
which may be enunciated thus.* 


Let B > «, f(s) be regular in the strip «a <o< f; 
if(s)| < C,(|t|-+1)* for o =a, 
if(s)| < Cy(\t}|+1)° for o= 8, 
if(s)| < Cx(\t}+1)% for «a<o<f, 
where C,, Cy, C3 are independent of s, a and } absolute constants. 
Then for a<o<fB 
If(s)| < 20, (|t|+-1) 


For the application, let 


a Ee, 
x= I+, 
=F = 3-4 
f(s) = Lis, x), X Xo 
In view of (20) and (21) and the easily-proved inequality 
| L(s, x)| < C,(\t}+1)% («a<o<f), 
it follows that the theorem applies, and therefore, on o = } 
\L(8, x)| < Ak*«(|t| 41). 
By using f(s) = k-**L(s, x), (22) can be replaced by 
|L(s, x)| <Akt+<(\t)-+1) 
For x= x9, remembering that x is primitive, and that therefore 


K=1, 
| L(8, x)| = |£(8)| < A(t) +1)** < Akt +<(|t]+-1)#+€. 


! 


3. Proof of (7) 
We now return to the right member of (17). On o= $+, 
ri. ee, 


[P(s)| << Ae"(|t|+-1)s; 


us, x)= [Ls x) Be) < ABE (|e +1), 
by (18) and (12); | 


\E(s, x,d)| < Ade < Ak. 


* Cf. Vorlesungen, ii, p. 49. 
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Accordingly, by (14), 
[ y-*T(s) > = Me” ds 


| 
| y-T(s)n(s, x)€(s, x.) ds}, (23) 


(4+e 
a0 


[ y~*T'(s)n(s, x)é(8, x, d) ds| < Aki +€|y|-#-€ [ e-taa—lare ylgt+e dt, (24) 
(t+e) 0 
Putting y = }7—|arc y|, then* 
0<y'<nly|, 
so that the right member of (24) becomes 


< Akt+¢|y| 


t-t»,-t-¢ e“udu 
0 
< Akt*ly|ttent+e, 
Substituting into (23) and making use of} 
| B(x)| < Akt+s, 
we get finally 


a(m)p™ 


ds < A(kly|)ient+ 


< Anité€, 
thus proving (7). 
4. The principal result 
Starting with (7) we prove 


Vyv-1 — 
= DI S(n) < Anv-t+e, 
as : 


where S(n) = > ta > am, 


k=1 pik) 
Since the proof is almost exactly like that of Theorem 249 of Vor- 


R,(n)— 


lesungen, i, it is unnecessary to reproduce it here. 
It remains, then, to identify S(n) with S(n). To begin with, it is 
easily seen that, if we put 


D(k) — 2 > om: 


plk) 
D(kk2) = D(k,)D(ke) for (k,, ks) = 1. 


* Cf. Vorlesungen, i, p. 213. 
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This enables us to write 


S(n) = 1 1+ D@)+D(p*)+ 


-_ {3 gies - 
| 0 for e>1 awe 


p’—p* for e<f 
—p* for e=f+1)} (p!|n, p+ / n). 
0 for e>f+1 
Substituting into (26), we see immediately that 
S(n) = S(n); 
this completes the proof of the 
THEOREM. For v >3; R,(n), S(n), b as defined in the Introduction ; 
A depending only on «, « > 0, 
bn” 


Bm) — | yy! 1 5) < An’-t+e, 


5. An application to the representation of a large integer 
as the sum of v ‘quadratfreie’ numberst 


We suppose that a(n) = 2(n). 


Then by (4) 


one 
pt)’ 


p* 


ani . “a 9. ani euteiit Zeal _ 6 
7(p*) = 0 for e> 2; b= [Te sl = 


S(n) can now be easily evaluated: 


[10 -T1 (yay) 
[, -T1(- dees) 


pin, p* in 
* Vorlesungen, i, Th. 220. 
+ Cf. Linfoot and Evelyn, ‘On a Problem in the Additive Theory of Num- 
bers’: Math. Zeitschrift, 30 (1929), 433-48, for a treatment of this problem 
along the lines of the Hardy-Littlewood solution of the Waring problem. 
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II " l] (1+ atnay=3) (I or) 
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oy 

, ~—r / (—1) : — ( 2__])r-1 
S(n) = || (1— aaa | | ip 
(p?—1)” 


This expression is obviously never zero. It is interesting to notice 
that the finite product depends only on primes appearing to at least 


the second power in n. 
6. A second application 
We shall here take 
a(n) = oi) r>}, 
n” ‘: 
o,(n) denoting the sum of the rth powers of the divisors of n. Then 


—(e+))r+1_} —(e+2)r+2_) 


by (4) p-?-+-p 


7p!) = 1+p?-1+ p-*-2.... 
=p, 


Substituting into (3), 
p! 


S(n) : | | (a) | | (1+ 25 ++ s) 


p 
_% v1 (2) 1 
yirt+bv-1 C{(r+1)}- 

ne 


Since rae 1 . 
b=| I(} | gait seat) 








THE SUMMABILITY OF A FOURIER SERIES BY 
LOGARITHMIC MEANS 
By G. H. HARDY 
[Received 18 February 1931] 


’ 1. THERE is a simple necessary and sufficient condition for the Cesaro 
summability of the Fourier series of a bounded function. Suppose as 
usual that f(t) is integrable in (—z, 7) and that 


f(t) ~ day+ > (a, cos nt + 6, sin nt). (1.1) 
I 


Further, suppose that f(t) is bounded near ¢ = x, and that the question 
is that of the summability of the series, for {= 2, to sum s. Then 
the condition t 
| d(u) du = 0(t) 
0 
or g(t)=0(1) — (C,1), 
wheret>Oand  4(t) = d{f(x+t)+f(x—t)—2s}, 
is necessary for summability by Cesaro means of any order, and 
sufficient for summability by Cesaro means of every positive order.* 

There are simple bounded functions which do not satisfy this con- 
dition. Suppose, for example, that x = 0 and 

f(t) = cos(a log |t}), (1.3) 
where a >0. Then b,, = 0, while a, behaves like 
A cos(a log n)+ Bsin(a log n) 


> 


nv 


where A and B are constants, and > a, is not summable (C). The 
condition (1.2) is naturally not satisfied; all means of f(t), round ¢t = 0, 
oscillate, substantially like f(¢) itself. 

The series } a, is, however, summable by Riesz’s logarithmic 
moans;} i 8, = fay ta,+0,+...+4,, (1.4) 
then : 

Ss 8 8 “ 
= (a+$+.+2] (1.5) 


logn logn n 


tends to a limit s when n-—> 00, in which circumstances we write 


* Hardy and Littlewood (3: see also 4). 

+ For the general properties of Riesz’s ‘typical means’ see the tract by 
Hardy and Riesz, and various papers by Zygmund in the Math. Zeitschrift 
and the Bulletin de l Acad. Polonaise. 
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lant Da,=s (R,1). (1.6) 


These facts suggest that we should try to find some general theorem, 
of the type of that just quoted, concerning the Rieszian summability 
of the Fourier series of a bounded function. Actually Theorem A 
below covers a rather wider class of functions. 

Some theorems concerning the Rieszian summability of. trigono- 
metrical series are known already. In particular Zygmund* has 
shown that the Fourier series of any integrable function, bounded 
or not, is summable (R, 1) whenever (1.2) is satisfied; and he and 
Jacobt have investigated more general theorems of the same type. 
But none of these theorems covers such a function as (1.3), or is 
quite of the type desired. 


2. THEOREM A. Suppose that 


t 
O(t) = | \d(u)| du = o(¢ log 7) 
J 


(a condition satisfied whenever (t) = 0 (log -) and in particular when 
f t 


f(t) is bounded near t= x). Then a necessary and sufficient condition 
that the series should be summable (R, 1), for t= x, to sum s, is that 


_ fem, fod 
w= | - lu - o(Iog 7) 
t 


when t>0. 

It is easily verified that (2.2) is satisfied whenever (1.2) is satisfied, 
and that the function (1.3) satisfies (2.2) but not (1.2). In proving 
the theorem we may make the usual simplifications, supposing that 
f(t) is even and that a, = 0, x= 0, s= 0, so that 4(t) = f(é). 


* Zygmund (7). 

+ Zygmund (8), Jacob (5). 

t Though in other respects they are more general, since they apply to series 
which are derived series of Fourier series but not Fourier series themselves. 

It may be useful to recall that (1.2) is not, when f(t) is otherwise unrestricted, 
a sufficient condition for summability (C, 1), though it is sufficient for sum- 
mability (C, 1+-8), for any positive 4, or (in virtue of Zygmund’s theorem) for 
summability (R, 1). A series summable (C, 1) is necessarily summable (R, 1), 
but a series summable (C, 1++-8), for every positive 8, is not, so that Zygmund’s 
theorem is not a corollary of the known results concerning Cesaro summability. 
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(i) Proof that the condition is sufficient. It is known* that if (as 
here) f(x) is even and integrable, and a, = 0, then 


x(x) =  f(u) 4 cot 3u du 


is also even and integrable, and 
x(x) ~ 4bo+ > b, cos ne, 


_ ay+a,+...+a,_,+4a, 
n 


where = 


for n >0. It is plain in the present case, from (2.2), that 


x(x) = o(log ) 


for small positive x. 
We write 
8, = 4,+4,+... Fy; O,, = 8, +8s+.--+8n-1 
1 — _— | | ‘ ‘ 
to = 4bo, ty, = 4b9+0,4+...+6, (n>), 7, = ty +t, +...+4,-4. (2. 
Then 


sin?3nx r sin*dnx 
rsaawell need g) 2" dx. (2.43 

sin?3x *) sin*3x ) 
From (2:32) and (1.5) 


Sn _ b S,,=t, +0 (log n). 
n 
From (2.33) and (2.43), 


1/n 7 
on [ o(1og=)n* dx re. [ (toe) S 
7 x w . x) x 


oy 


0 1/n 


7) (m = log n| +o(nlogn) = o(nlogn). 
n 


Combining this with (2.42) and (2.5), we obtain 
S,+8,+...+8, =7,+0(nlogn) = 0(nlogn). (2.6) 
But (n+1)8,—S,—Sg...—S,, = 8 +84... +8n; 
and so, by (2.6), 
S, = as ws Cn 


n+] 


oe =a 1 t0(ogn). (2.7) 





* Hardy (2). 
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On the other hand 


” So ae 
Sn 4 [payne de +0(1) = J+0(1), 
n x 

0 


say, and 


ned fi 


n 


sin? }na 
x) | dx 


say. Here 1/n 


| f(x)| dx = "(=)=0 o (log n), 


a 
« 


0 
by (2.1), and - - 
& (x : P(x) dx-+-o(1) 


>= = — ; 
na x nar J x 
l/n 


= o(log 7). 


Collecting our results, we see that 
o (log n), 


and so, from (2.7), S,, =o(logn), the result required. 
3. (ii) Proof that the condition is necessary. We now assume that 
, = 0(logn) and deduce (2.2). We begin by showing, without using 
2.1), that 1 
x(#) = 0 (log -) (C, 2), (3.1) 


i.e. that ¥.(z) = du x(v) dv o(2*log ). (3.11) 
x 


The proof is much like that of Riemann’s classical theorem. We 


a - a 
9 X2(*) 1, LSy sin }na\? 
a 2 - 9 IoT dn 1 = 
”, “ — ane , 


x 
1 


have 


= t,.{g(nx)—g| (n+ 1)2x]}, 
sin *) ; 


where g(u) = 
i. 





THE SUMMABILITY OF A FOURIER SERIES 111 


and so 
(n+1)z - 


2 tage) <> Il J ig’ (u)| du => |t,\I,, (3.2) 
Ms 0 
say. But g’(w) is O(u) for small w, and O(u-*) for large uw; and so 
I, = O(nx*) (nx < 1), i,= o( : ) (nx > 1). 
nx, 
Substituting in (3.2) and observing that, after (2.5), t, = o(logn), we 
obtain 


> o (log n)O(na?) + ba o (log n)o( : 
nx<l nx>I1 “als 


=: o(2* 2 n log n)- +o(~ Z a) o(lox. ): 


nw s 


which is (3.1). 
It follows that (x) = 0 (2° log) (3.3) 


(suffixes denoting, as before, integrations from 0 to x). On the other 
hand it is easily verified* that 
Je%p(x) = p(x) —4af,(x)—3fa(2). (3.4) 
The last two terms are o(x*log-), by (2.1); and it therefore follows 
x 


from (3.3) that ' 
Hx) = o(log =) 
which is (2.2). This completes the proof of the theorem. 

4. We might naturally express (2.2) by saying that f(t) is ‘con- 
tinuous (R, 1)’ for t=2; Theorem A then asserts that, if (2.1) is 
satisfied (and in wiationiion if f is bounded near t= x), continuity 
(R, 1) ts a necessary and sufficient condition for summability (R, 1). 

It is natural to ask what happens when the condition (2.1) is 
dropped, and the answer suggested by the analogy with Cesaro 
summability is as follows. We must begin by defining summability 
(R, k) and continuity (R, 1). We may then expect a double scale of 
theorems like those investigated first by Littlewood and myselff, and 
then, more precisely, by Bosanquet and Paley.{ None of these 
theorems, however, will be of the ‘necessary and sufficient’ type; for 

* For example by differentiation. + Hardy and Littlewood (3). 

t Bosanquet (1), Paley (6). 
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that, it will be essential to close the cycle by some such condition as 
(2.1). It might be worth while to push the analysis a little farther 
(though hardly to develop it in full detail), but I confine myself to 


what seems to be the simplest and most interesting case. 
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By MARY L. CARTWRIGHT 
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1.1. IN a previous paper} I considered functions of the form 
1 
f(z) =fle+iy) =fire!®) = [ e4(t) dt, (1.11) 
“1 
where ¢(¢) is a complex function, integrable in the sense of Lebesgue, 
and ¢(t) tends to a finite limit other than 0 at each end. The object 
of this paper is to consider the cases in which ¢(¢) tends to 0 or 00 at 
one or both ends. If 4(t)-> 00 as t> +1, we need lighter restrictions 
in order to obtain equivalent results, and if 4(¢)—> 0, heavier ones, as 
the following example shows. 
Ss se the 
uppose that b(t) = 1—t (0<t<}) 
=I14+# (—1<t<0), 
t ? 
° 4sinh?(3z 
then d¢(t) = | ¢'(u)du, f(z)= = ) 
As 1 
Hence, although ¢(/) is an integral, the zeros of f(z) occur in pairs, 
and not as in Theorem VI of Z. 
In the above example all the zeros still lie on the imaginary axis, 
but Titchmarshf{ has constructed a function of the form (1.11) which 


29 


~ 


has an infinity of real zeros. 
1.2. I suppose first that 
d(t) ~ C(1—t)4 as t>1, 
d(t)~C’'(1+t) ast>—1, 
where g> —1, C40, C’ +0. We may suppose without loss of 
generality that C = C’ = 2%, so that we may write 
p(t) ~ (l1—#) (1.21) 
as t+-+1. The other cases in which C 40, C’ + 0 may be reduced 
to this by simple transformations as in Z, § 1.4. 
Theorem I gives asymptotic values for f(z) in the regions 
lcos | >8> 0, r > 1r9(5), 
+ M. L. Cartwright, Quart. J. of Math. (Oxford), 1 (1930), 38-59. I shall 
refer to this paper in future as Z. 
t E.C. Titchmarsh, Proc. London Math. Soc. (2) 25 (1926), 291, Theorem ITI. 


3695.2 I 
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and shows that f(z) has no zeros in those regions. This is the only 
theorem deduced from (1.21) alone. We seem unable to improve on 
the formula dish Sain, 
which was proved by Titchmarsh? for any function of the form (1.11), 
without further conditions on ¢(t). 

When gq is an integer, we impose the conditions used in Z on the 
qth differential coefficient of ¢(¢), instead of on ¢(¢) itself. If gq is 
fractional we have a complication arising from the irregularity of the 
dominant term itself; for if g = m-+a—1, where m is 0 or an integer 
(0 <a < 1) then the mth differential coefficient of 4(¢) becomes infinite 
at ¢—-+1. Theorems V and VII give typical simple cases for g = a—1; 
the cases in which m>0 reduce to this on integrating m times. 
Theorems VIII and IX are the most general results of their kind; 
they include Theorems II, III, V, and VII. Generalizations of IV 
are also indicated briefly; I and VI stand by themselves. The use of 
integrated Lipschitz conditions gives very great generality, and is 
equally suitable whether q is an integer or not. By separating the 
term (1—/?*)? from the rest we avoid confusing the well-defined 
irregularity of (l1—?t)*-1, where 0<«< 1, with the general irregu- 
larity of d(t). 

In §§ 7.1—7.3 I use the methods evolved in the previous theorems 
to prove a new type of theorem. Returning to the case g = 0, con- 
sidered in Z, I suppose that ¢(t) belongs to Lip(«, p), where ap > 1. 
The error term for n(p) in these results is O(r°-~»’), where 


p’ = p/(p—1), 
while a,, the nth Fourier coefficient of ¢(t), is of the form O(n-*). 
Hence the analogue observed in Z{ fails here. Some results can be 
deduced directly from the Fourier coefficients, but only if a,, is o(m-*) 
or even smaller; these, however, I must leave for another time. 

In all these cases the zeros are distributed more or less like those 
of 2-7-1 J,,, 
yz. 

These methods can also be applied when 


(—iz), ie. like those of z~¢—sinh(z-+-}q7i) according as 


p(t) ~ C(1—t)4 as t>1, 
p(t) ~ C'(1—t)¢ as t> —1, 


+ E. C. Titchmarsh, loc. cit., Theorem IV. t See § 3.1. 











THE ZEROS OF CERTAIN INTEGRAL FUNCTIONS 115 
where g ~q’. We reduce as before to the form 
P(t) ~ (1—t)(1+4)". (1.22) 
The zeros of the typical function are given approximately by 
Nqg+1) , 


a ae 9 Si te 
x= h(q—q' )log nat blog oti i) “a (1.23) 


y= +[r+hq+q) +4, 
where e¢,, is a number (not necessarily the same in both places) which 
tends to 0 as noo. When ¢(t) is of the form (1.22), the zeros of 
f(z) approximate to (1.23); and owing to the lop-sided arrangement 
of the zeros with regard to the imaginary axis we have to impose 
heavier conditions on ¢(¢) at one end than are needed at the other. 
This makes it difficult to state satisfactory theorems. I shall there- 
fore only give a brief sketch of the chief points of difference between 
these functions and the previous ones. 
Finally, I discuss briefly functions for which ¢(¢) tends to 0 ex- 
ponentially at the ends of the interval of integration. 
1.3. I use the notation of Z as far as possible, that is to say: 


z, = 1,e%, z= ree™,..., where O0< 1, <1. <..., are the zeros of f(z) 


n(p) is the number of zeros for which |z| < p; N(p) = [ a du; a,, is 
u 


« 


0 
the nth Fourier coefficient of 4(¢); [p] denotes the integral part of p. 
[ am indebted to Professor Hardy for many valuable suggestions 
particularly with regard to §§ 4.1—4.4 and § 9.1. 
THEOREM I. Jf d(t) satisfies (1.21), then, given any « > 0,85 > 0, 
we can choose 14 = 19(5, «) 80 that 


~4 +1 
1+1f(z) — e-2-a7i( + £) (—39r#t+8 <0< —47—58) 
201(q+1)— : : 


= (140) (—Jn48 <0 < 4-8) 
=es7i(14f) (fn +8 <0 < fn), 
where € is not necessarily the same in each case and || <e« for every 
r >, 
CoroLLtaRy. If d(¢t) satisfies (1.21), then f(z) has no zeros in the 
regions |\0| < $a—8, |r—O0| < 4n—8, r > 19(5). 
It is convenient (but not necessary) to use these asymptotic expres- 
sions to prove some of the later theorems. The apparent difference 
12 
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between the first and third is only due to the many-valued factor 24 
on the left-hand side. 

Proof. Suppose first that x > 0. We may write 4(¢) = (1—#?)"p(t), 
where (t) is continuous and equal to 1 at ¢=1, and, given any 
€, > 0, we can choose 7 so that 


ls, (t)| = |22—(1+4)%f(1)| << «, (2.11) 
for l1—y <t<1. Hence 
1 J 1-9 
fe) = 24 | et(1—tyrdt + | e(1—t)r y(t) dt + | eg(t) dt 
i—y i—y —1 
= J,+1,4+-Js. 


It is easy to see that 
hy 
I,| < e-e | p(t)| dt< Ae-v2,t 
—1 


1 7 





while I,| <<, ( e(1—t)4 dt = e,e* [ e~ta dt 
1—n 0 
Py or] ' = 1 
< ee? [ emerdt =e, GT) (2.12 
3 gat 
0 
” : . a 
Also I, = 2%e | e-*0 dt aves) | — [ jesen at 
0 \; 7 


- 24e°T (q+ 1)2-4-14 Ofer}. 
Putting these results together, we have 

2ittf(z) = (q+ 1)2%*(1+-Z), 
where |C| < «,(cos@)-4-1!+ Art+le-™ < e, provided that |6| << jn—8, 
r>r(e,5). For we can first fix 6, and then choose ¢, so that 
e,(sind)-¢-1 < $e. We then choose » so that (2.11) is satisfied, and 
lastly 7) so that Arg@He-10 sind — Le, 

The other two asymptotic values are obtained in a similar fashion 


by dividing the range of integration near t = —1. 
3.1. Suppose that g = m, where m is a positive integer. In the first 
example of § 1.1, d(t) ~ 1 t as t> +1, and ¢’(t) is of bounded varia- 


tion, while the zeros are like those of the example in Z, § 3.11, in 
which ¢(¢t)—1 at the ends and 4(¢) itself is of bounded variation. 
+ Compare the formulae for Bessel functions, see G. N. Watson, Theory of 
Bessel Functions, Cambridge (1922), 48. 
{t Throughout the paper I use A to denote a positive constant, not neces- 
sarily the same in different places. 
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This suggests that we must now impose conditions on ¢™(t) instead 
of d(¢). When this is done, we obtain results similar to those of Z, and 
the example shows that we cannot expect anything much better. 
The continuity condition at the ends corresponding to that used 
in Z is that ¢™(t) is continuous. Suppose that this is satisfied; it 
then follows from Taylor’s theorem that 


f(t) = o(1)+(¢—1)¢’(1)+...4+ omen + €(t—1)}, (3.11) 
where 0 < €< 1, neart=1. ciunasiadiies ies and (3.11) we have 
$(1) = $'(1)=...=¢™-(1) =0, (3.12) 
and ¢(1) = m!(—2)". Similarly 
$(—1) = $(—1) =... = ¢"-(—1) (3.13) 
and 6™(—1) = m!2™. 
THEOREM II. Suppose that ¢(t) is of the form (1.21), and that 


f(t) is an integral. Then the zeros of f(z) are determined asymptotically 
by the formula 


z= +(n+4m)zi+e,, (3.21) 


where n is a positive integer and «, 0; also 


n(p) = 2|2—am|, (3.22) 


7 
provided that p is large and that p/m—4m is not too near an integer. 
We can integrate m-+-1 times, and we have 
1 
( ikea " ‘ 
1 ° 2t (m)(¢ C 
J *) zm etd ) dt 


off 


1 
ais ; — m ff 3 
a ” {e* gm™(1 —é “g(— 1 i = 4 eign M(t) dt 
v = 9m 
: ; Le {e*— ( —ipe-4—f. 
But, by Titchmarsh’s lemma,f J = 0 {e*'/r™*4}, so that 
gm If(z (z 
D(m+ - 


where «, denotes generally a function of r and @ which tends to 0 


e*(1+e,)—(—1)™e-*(1+-«,), 


uniformly in 6 when r—> oo. 
At a zero e2 — etmmi(] te), 
so that 2: = i ; 


+ E. C. Titchmarsh, loc. cit., L nma 2.2, or Z, Lemma A. 
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We can complete the proof of (3.21) as in Z by showing that one 
and only zero is associated with each large integral value of n. We 
prove (3.22) by showing that z”*1/(z) has the same number of zeros 
as sinh(z—4mz7) inside |z| = p, provided that the circle does not pass 
too near a zero of sinh(z—4mz?t). 

3.3. THrorem III. Suppose that 4(t) is of the form (1.21), and that 
h(t) is of bounded variation and continuous at t= +1. Then 

(i) all the zeros of f(z) lie in the strip \x|< K for some K; 

9 
(ii) n(p) = + O(1); 


7 
9 
(iii) N(p) = =o — (m+ 1)log p — log|f(0)|-+-0(1). 
We integrate m times and then apply the methods of Z, §§ 3.12— 
3.42.¢ The results follow in precisely the same way. 
3.4. Similarly we have 
THEOREM IV. Suppose that 4(t) is of the form (1.21), and that d™(t) 
is continuous with modulus of continuity w.t Then 
(i) all the zeros lie in the region \x| < Krw(r-), 


(ii) m(p) =P + Ofpe(p~)}, 


(iii) N(p) = “P — (m+ 1)log p — log |f(0) |+ Ofpw*(p-1)}+- 


+0lw a || a \. 
\ (p ee) 


p(t) ~ (1—#)*"*, (4.11) 
where 0< a< 1, as ¢-> +1, and consider first the case corresponding 
to Theorem II. A simple transformation shows that }f(nzi) is 
actually the complex Fourier coefficient of ¢(¢/7). Bromwich,§ 
Young,|| and Haslam-Jones** have obtained asymptotic values for 
the Fourier coefficients of different classes of functions satisfying 


1 
+ On p. 43 of Z in the last line, read f |\dp| insteadof max  ([4(t,)—¢(t,)]- 
1-8 1-8<thi<te<1 


t Ch.-J. de la Vallée Poussin, Legons sur l’approximation des fonctions d’une 
variable réelle, or Z, § 4.1. 
§ T. J. ’a Bromwich, Theory of Infinite Series (2nd ed.), 494. 
W. H. Young, Proc. Royal Soc. (A), 93 (1917), 49-55. 
** U.S. Haslam-Jones, Journal London Math. Soc. 2 (1927), 151-4; see the 


4.1. Now suppose that 


note added at the end. 
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(4.11). In each case the working holds equally well for non-integral 
n so that we can easily obtain from it an asymptotic value for f(iy). 
From this we deduce the required result by a theorem of Phragmén 
and Lindeléf. This method is not directly applicable when g = 0 or 
an integer; for we have to go through the working in order to find 
out the dominant terms. For example, if g = 0 and ¢(¢) is an integral, 
pin~ 24 0(2), 
vy Y\, 
although, since e”7' = e-"*', we have a, = o(n-!). A mere order con- 
dition on f(iy) is insufficient here; we must be able to separate the 
principal terms from the rest. 
4.2. Using Bromwich’s work, which is the simplest of the three, 
we have 
THEOREM V. Suppose that $(t) = (1—é)*—"f(t), where p(t) is of 
bounded variation in (—1, be he (l—»,1), and p(t)>1 as 
t>-+1. Suppose further that 4(t) belongs to Lip* at in (—1+-y, 1—7). 
Then the zeros of f(z) are determined by 
z= +{n+}(a—1)}vi-+e,, (4.21) 
where n is a sali ta integer and «, >0. 
We have 


fly) = 


1 
+f] \ cith(t) dt = Iy+Iy+Ty. 
—1 -1liyn 1- ,! 
Since (1+-4)*-(é) is of bounded variation in (1— », 1), it follows from 
Bromwich’s work that 

7 

I, = e'¥ | e~‘ugo-1(2—t) 9-1 —t) dt 
0 


= DP(ae'v2*-Miy)-*+-0 (|y|-*). 


Similarly I, aed P(a)e-¥2°-1(iy) “Ae t+amt to ( y x). 


Also , ” 1—n+mly 


. 


i= | e'v'd(t) dt = — [ eimg(t— “) dt 


/ 


ln —l+ntaly 

1—9~aly s 7 ” 
} | eivt hag -4(t— =) dt +43 @ ie iP e'V'd(t) 
Lintaly 1-- a 


+ d(t) belongs to Lip « in (a, b) if eee = O(h®) (a <t<tth<b), 
and to Lip* « if |f(¢+h)—d(t)| = o(h®). 
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Hence 


” 
\Te| <3 { P(t )-4(*—2) dt 4 
a 7 
= o(ly\-*)-+O(\y|-*) = o(\yi-). 
We now have 


F(z) = z%e-*f(z) —T'(a)2%-1(1 —e-2+ 27%) — 9 (1) (4.22) 


for @ = 37, and also, by Theorem I, for 6= 0. It is easy to see that 
F(z) = O(r%) uniformly for 0 < @ < 37, and it follows from a theorem 
of Phra; ugmén and Lindeléf ft that (4.22) holds uniformly for 


0<@< 4a. Using a similar process for }1 < 0 < 7z, we have 


fie) = Te ee  ¢-+171) 40%) 


ri ore 4 
Na 7 


for y >0. Hence at a zero in the upper half-plane 
e*(1 re)= = ee Ha—Dri(] re 
so that z= {n+}(a—1)}n 
In the lower salt aiid we have 
Pee, ees ee 
f(z) = = te spices )+0(—), 
so that z= —{n+}(a—1)}i+e,. 
4.3. In general we have 
THEOREM VI. Suppose that $(t) satisfies (4.11) with O0<a<1, and 
that (iy) *f (ty) = T"(«)2% lielv | e—tys amt) 9(1), (4.31) 


for y20. Then the zeros of f(z) are given by (4.21). 

It is easy to show that functions for which 4(¢) satisfies Young’ st 
or Haslam-Jones’s§ conditions at = +1 satisfy (4.31), and there- 
fore have their zeros given by (4.21). Young’s condition at the point 
t =a is that 


(t—a)?-%'(t) is bounded in some interval (a—n,a+) and tends to 
las t—a. (4.32) 
Haslam-Jones’s condition is that 
h(t) = (t—a)*—b(t), where y(t) > 1 as ta, and y(t) satisfies ‘ Young’s 
a+t 
convergence criterion’ ,|| i.e. [ |d{(u—a)b(u)}| = O(t). (4.33) 
a 
+ G. Pélya and G. Szegé, Aufgaben und Lehrsitze, » vol. i, Berlin (1925), 149. 
ee se cit. § loc. cit. 
. W. Hobson, Theor y of Functions of a Real J lesa, 3rd ed., 2, 532. 
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4.4. Corresponding to Theorem III we have 


THEOREM VII. Suppose that $(t) satisfies the conditions of Theorem 
V except that it belongs to Lip « instead of Lip* « in (—1+ 7, 1—7). 
Then (i) and (ii) of Theorem III are satisfied, and 


9 
(iii) N(p) = © — alogp + log|f(0)|-+0(1). 
7 
We write 1 ~1+9 1—» 
}e“g(e dt + | e“d(t) dt 
~1 —1+7 
= f(z) +fo(2). 
As in Theorem V, we have 
7 9a-1 x 
* I (x)2 a (e?+-e “~+0(—) (4.41) 
= f= 
for y 
1—n+7/y : 
[ eft4 ivMd(t) dt — | eft + iy ering (t—") dt 
—1+ 1 tarly , 
i—y 


J 


—1+n+7/y 


t= | a+ 
7 


efx iyt 40 — “ening 
14 


( Bi, 1—7 


+] h(t) dt 
iJ me | 


= J,+/,. 
It is easy to see that I,| << Ae?-2/y. (4.42) 
Also 


I,| <e-v2 f $)—$(t—) +(e emuy(t—2) at (4.43) 
* 7 


1—7 7 1—7 
ef [ 4(t) -$(t—7) dt +el-m2(1¢-27) [ \b(t)| at 
ee Y fea 
A / el—nx 1A ae —n)a ie (4,44 ” e(l—nix 
y* i y* 
as yoo. Hence we can choose K so that 
\z"e-*f,(z (A,+Agr)e-™ < } (4.44) 
for «= K. It follows, as be a from the theorem of Phragmén and 
Lindelof that (4.44) holds wen for x > K. We now have 


zf(z) ee- 1¢7(1+-Z), 
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where |¢|< 1 for x > K. Then (i) follows immediately, and it follows 
from (4.41), (4.42), and (4.43) that f(z) = O(r-*) for all fixed x. We 
can then apply the methods of Z to z%f(z), and obtain (ii) and (iii). 
4.5. If d(t) = (1—#*®)*-4h(6), where A(t) is continuous with modulus 
of continuity w, we have results corresponding to those of Theorem 
IV. The proof is like those of Theorems X, XI, and XII of Z, except 
that we write 4(¢) = (1—t)*-12%-1+-4,(t) in the end interval. 
4.6. Theorems V and VII are by no means the most general of 
their kind, and the abrupt change in the hypotheses at the points 
=-+(l—y) is quite artificial However they are included in 
Theorems VIII and IX. Theorem VI, although very general, is 
unsatisfactory without some general criterion as to whether the 
hypothesis is satisfied for any given ¢. 
5.1. Suppose that 4(¢) is of the form (1.21 


m is zero or a positive integer, and 0< «<1. Suppose also that 
qm 
dti™ 


where now g = m+a—l, 


gmt) = (1-2) 0(1-2)9, (5.11) 
ast+>-+1. This is equivalent to the continuity conditions used in 
§§ 3.1 and 4.1, and it follows that (3.12) and (3.13) still hold for 
0<a<l. 

Consider first the case corresponding to Theorems III and VII. 
Hardy and Littlewood} have shown that functions of bounded varia- 
tion are equivalent to functions which belong to Lip (1,1). It there- 
fore seems natural, when a < 1, to replace the hypothesis of bounded 
variation by ‘d(¢) belongs to Lip («, 1)’.+ The definition of Lip («, 1) 
is as follows:t Suppose that 4(t) is integrable and periodic with period 
b—a, and suppose that 


+h)—d¢(t)| dt = O(\h|*), (5.12) 


where 0< a<1, as h+0; then d(t) belongs to Lip («,1) in (a,b). We 

say that 4(¢t) belongs to Lip* (a, 1), where 0 < a < 1, if the above defini- 

tion is satisfied with o instead of O. If it is satisfied for « = 1 with o, 

then there is a constant c such that 4(¢)=c almost everywhere.§ 

For Fourier coefficients, 4(¢) an integral gives an o-result compared 
. H. Hardy and J. E. Littlewood, Math. Zeit. 28 (1928), 619. 


. H. Hardy and J. E. Littlewood, loc. cit., 612. 
E. C. Titchmarsh, Journal London Math. Soc. 2 (1927), 36, 37. 
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with an O-result for ¢(¢) of bounded variation. For the purposes of 
this paper, therefore, I shall say that ¢(¢) belongs to Lip* (1,1) if it 
is an integral. We now have perfectly smooth generalizations of the 
properties of bounded variation and of being an integral which are 
applicable for 0< «<1. We say further that 4(¢) belongs to Lip (q, 1), 
where q = m+a—l1, if 


gm -2(t) = i g(t) dt, (5.13) 
ig 
and $™(t) belongs to Lip (a,1). We may of course replace Lip by 
Lip* in this definition. 

The relations of the integrated Lipschitz conditions to those used 
in §§ 4.1—4.4 are given by the following lemmas. 

LemMA A. (t) belongs to Lip* (a,1), where «< 1, in (—1,]1), if tt 
is defined outside the interval by periodicity and satisfies one or other 
of the following conditions inside the interval: 

(1) d(t) = (t—a)* Y(t) (—1 <a<1), where (t)>0 as t>a, and 
u(t) is of bounded variation. 

(2) (t—a)?-%'(t) is bounded and tends to 0 as t>a.t 

(3) d(t) = (t—a)*(t), where p(t)>0 as t>a, and y(t) satisfies 
Young’s convergence criterion.§ 


(4) A(t) = (t—a)P(t), where B> a, and x(t) belongs to 
Lip* («+ 1—f). 


(5) A(t) belongs to Lip* « 

The last class of functions obviously belongs to Lip* (a,1), and 
the proof for (4) is very straightforward. (3) includes (2) and (1), 
and the proof is a little more difficult. 

Lemma B. ¢(t) belongs to Lip («,1) in (—1,1), if it ts defined outside 
the interval by periodicity, and satisfies one or other of the following 
conditions inside the interval: 

(1) A(t) satisfies (1), (2), or (3) of Lemma B except that y(t) o 
(t—a)?-%d'(t) may tend to a limit other than 0. 

(2) h(t) = (t—a)Pg(t), where B > «, and y(t) belongs to Lip (x+1—B). 

(3) A(t) belongs to Lip «. 
5.2. It is easy to verify that, when a < 1, (1—?*)*~* does not belong 

+ Compare Theorem V. 


{ Compare (4.32). If #4) ) satisfies (4.32), d(t)—(t—a)*"! satisfies Lemma 
A (2). § Compare (4.33), and see footnote. 
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to Lip*(a,1), but only to Lip(a,1). Hence, when ¢™(¢) is of the 
form (5.11), it cannot, in general, belong to Lip* («,1). We therefore 
subtract (1—/*)* from ¢(¢) before imposing our hypotheses. This 
makes no difference if ¢ is an integer. 

THEOREM VIII. “ar that $(t) is of the form (1.21), and satisfies 
(5.11), and that $(t) = (1—t#?)4+- x(t), where x(t) belongs to Lip* (q+1, 1). 
Then the zeros of f(z) are given by 

z= +(n+4q)mi+e,, 
where n is a positive integer and «, > 0. 

5.3. In the theorems about functions of bounded variation we only 
assume that ¢(¢) or d(t) is continuous at t= +1, but from this we 
deduce that ¢(¢) has continuous variation at those points,7 i.e. 


1 
lim dd 


= 1 ; n 


or, what is practically the same thing, 


( d(t)+ 


lim lim 
n—0 h a F 
—7 


We now require a similar property for functions which belong to 
Lip (a«,1), where a< 1. The natural generalization is as follows: 
(t)>0 Lip (a,1) as t>a+0, if 

H(t) = of(t—a)3, 


where 0< a< 1, as t>a-+0, and 


lim lim 
70 h->0 ¥ 
a 


Also ¢(t)->0 Lip (q,1), where g = m+a—1, if 
g(a) = $'(a) =... = $2), 


and d"-V(t)>0 Lip (a, 1) as t>a+0. 

It might be supposed that if 4(¢) belongs to Lip («,1), then (5.31) 
implies (5.32) for a< 1 as well as for a= 1. But the proof for «= 1 
depends on the expression of a function of bounded variation as the 


+ See footnote to § 3.3. 
t Compare G. H. Hardy and J. E. Littlewood, Math. Zeit. 28 (1928), 619. 
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difference of two monotonic functions. There seems to be no corre- 
sponding result for functions which belong to Lip («,1), where «< 1. 

Lemma C. ¢(t)—>0 Lip (a,1) as t>a+0 if 

either, (1) (t) belongs to Lip* (a,1) and satisfies (5.31), 

or, (2) A(t) satisfies (2) or (3) of Lemma B. 

THEOREM IX. Suppose that d(t) is of the form (1.21) and satisfies 
(5.11). Suppose further that d(t) belongs to Lip(q+1,1), and that 
A(t) = (1—#*)*+ x(t), where x(t)>0 Lip(q+1,1) as t>+1. Then (i) 
and (ii) of Theorem III are satisfied, and 


(ii) N(p) =P —(¢-+1)log p — log f(0)|+0 (1). 


Again the subtraction of (1—?é*?)* makes no difference if g is 0 or 
an integer, while for non-integral g, (1—#*)* itself does not satisfy 
(5.31) or (5.32). 

5.4. The dominant terms in the asymptotic value for f(z) are given 
by the following lemma: 


LemMA D. As r>oo 
1 
* a 24 ‘ et 
[eq —erye dt — PIED ee ¢-ssary 4 Qf | 
4 | ra+2] 
for 0<@0<z, while for —7<0<0 we have to read e-*-97 in the 


i +1 ai A 


second term. 
This is derived by easy transformations from well-known formulae 
for Bessel functions,} since we have 
8,(z) = J,.4(—t)P(q+1)P(3)(— diz) +. 
5.5. Proof of Theorem VIII. We have 
; 
flz) =s,(z)-+ | etx(t) dt (5.51) 
-1 
= 8,(2)+1(z). 
We define y(t) outside the interval (—1, 1) by periodicity. Suppose 
that ¢g = m-+a—1; the case a= 1 has already been dealt with. For 
x< 1 we have 1 
(—3P f 5 
T(z) =-—— | e@y(6) dt, (5.52) 
-1 


+ G. N. Watson, Theory of Bessel Functions, Cambridge (1922), 48. 
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where x(t) belongs to Lip* (a,1). Hence 
1 


1 
(iy)™I (iy) ae (—1)" [ eiut xm ) dt - =(— ])m+1 | emi (t+) dt 
i A y 


ont -1 


1 a. 
a} It. 


/ 


nme m . 
== ( [ emt xn == x(t + 
a 
It follows that 


a) eee 

I (iy) lim J X(t) —xX™|_ tH 4 " 
= of iy|-"-9} = of ly. 
Also it is easy to see that J(z) = O(e'*') uniformly in 6. Hence applying 
the theorem of Phragmén and Lindeléf to z%*4e-*J(z) in the right 
half-plane, and to z¢+1e*J(z) in the left half-plane, we have 
I(z) = 0 (e*'r-4-1) 
uniformly in @. 
It now follows from Lemma D that 

teIf(z) = D(g+1)24(e*—e-**4™) + 0 (e") 
for y> 0, and a similar formula holds for y< 0. From these we 
determine the zeros as usual. 

5.6. To prove Theorem IX we proceed as above. We deal with 
the interval (1—8, 1) by means of (5.31) and more especially (5.32). 
The rest of the interval is treated more or less as in § 4.4. 

6. We can easily generalize Theorem IV by using a function (8) 
corresponding exactly to the modulus of continuity. Suppose that 
w(5) 40 (5!-F), where 0<B< 1. Then ¢(é) is said to be continuous 
(8B, w), if |h(t.)—d(t,)| < 88a(8), provided that |t,—t,|)<6. If 
«w(5) = 0 (8!-F), we use the condition |¢’(t,)—¢’ (t,)| < 88 —w(8) instead. 
¢(t) is said to be continuous (m-+-f, w) if g(t) is continuous (f, w). 
Now, if ¢(¢) satisfies (1.21) and is continuous (q, w), we have results 
like those of Theorem IV. The methods are much the same as those 
used in Z. 

7.1. We shall now consider functions which belong to Lip (a, p),7 
i.e. functions for which 


b = 
[ |P(t-+h)—(t)|? ar| = Ofh}. (7.11) 


t G. H. Hardy and J. E. Littlewood, Math. Zeit. 28 (1928), 612. 
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If d(¢) belongs to Lip(a,p), where ap> 1, then it belongs also to 
Lip («a—1/p).¢ It follows at once from Theorem XI of Z that when 
(> Last >Th, nlp) = 2p/n-+O(p-+4P)., (7.12) 
We can, however, obtain a better result. 

THEOREM XI. Suppose that ¢(t) 1s of the form (1.21) with q=0, 
and that 4(t) belongs to Lip (a, p), where ap> 1, p> 1; then 

(i) all the zeros of f(z) lie in the region |x| << Kr“-™»’, where 


p’ = p/(p—}), 
9 
Oa 


If «= 1 and p> 1, ¢(¢) is an integral;t and therefore (5.21) holds. 
If poo, p’>1, and the theorem reduces to one for an ordinary 
Lipschitz condition.§ Since 
ap—l 
p(p—1)’ 
we see that Theorem XI (ii) is a better result than (7.12). 


7.2. Proof of Theorem XI. We can write 


1 
l1—a+——(1l—a)p’ = 
Pp 


sinh z 


1 1 
f(z) = [tay fexma—* —n oe E, 
“1 -1 


~ 


where y(t) belongs to Lip (a, p). We first find the value of J on the 
boundary of the region defined in (i), i.e. on 


f= Ky’, (7.21) 
The case «= 1 has already been proved in Theorem VIII, so that 
we only have to consider a< 1. In this case x>0o with r, and also 
y/%—> oo with r on (7.21). We shall make use of these points in the 
proof. 
Suppose that a> 0, y> 0; then on (7.21) we have 
1 


[I(a-biy)| < fe x) ~~ ase [+ {Jao 


y 


= I,+I-. 


. H. Hardy and J. E. Littlewood, loc. cit. 627. 
. H. Hardy and J. E. Littlewood, loc. cit. 27 (1927), 599. 
heorems X, XI, and XII. 
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Since x(t) belongs to Lip (a, p), it belongs to Lip (a—1/p); also y(t) > 
as t>+1, so that x(t) = Of(1—#2) 1p} (7.22) 
as t+>1. Hence 1 


I,< Ae | (1—t)*-? dt << — con 


ym" 1—-l/p 


Again 


4 dt 


-1 
= I,+/,+/;. 

Using (7.22) we can showft that J, < Ae*/(yx*-"”); and it is easy to 

see that [,< Axe™-”/y. Lastly, by Hélder’s inequality, 


! = at| a 
-y 


j 1 
<A f pp'act D ‘(J ? 
I, - | | é€ tf ) x(t | ~ glp’y : 


~1 

where p’ = p/(p—1). ¢ 2. these results we have 
Ae Arc A _* 

_ alp'yx | y Vga py ~ r 
for every r> 7, on (7.21), provided that K is sufficiently large. 
Applying the theorem of Phragmén and Lindeléf to ze-*J(z) in the 
region x > Kr@-~»’, we find that (7.23) holds throughout it. Hence 


I(a+ty)| - (7.23) ° 


y 


af(z) = e*(1+-¢), 
where |¢|< 1 for x > Kr@-~»’, From this (i) follows at once and a 
similar expression for « < —Kr“-~»"; and (ii) follows as usual. 

8. When d(t) ~ (1—t)4(1+2)”, where q > q’, the typical function is 

1 

84, ¢(@) = | e(1—t( +0" dt, 

—1 
and it is easy to show that 
C(q+ dal 


24+ 


P(g’ +1)2%-* 


aq +1 


$a) 1+e,)—(—1)" (1+¢,). 


The method used by Hardyt then shows that the zeros of s, ,,(z) are 
given by (1.23). 

Near t= 1 we want hypotheses suitable when ¢(¢) ~ (1—#)? and 
near t= —1 those suitable when ¢(¢) ~ (1—?*)”. In the less precise 


+ Compare (2.12). 
+ G. H. Hardy, Proc. London Math. Soc. (2) 2 (1904), 423. 
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theorems the middle of the interval may be neglected;} otherwise we 
must assume enough to make the hypotheses join properly. 


9.1. Finally we may ask what happens when ¢(¢)—> 0 very rapidly. 
There is one function for which calculation is easy. 
THEOREM XII. If o(t) = (1—t)-*?e-"4~, then the zeros of f(z) are 


quwven by t= (in ar) 2(1+-e,,) | (9 11) 
+na+(3nm)!(1+e,).J | 
We have 


1 2 
ean % _o f .# yp dt _ e-a- yy at 
(1 —t)3i? a= ae 
0 


« 


I,+-J,. 
It is well knownt that J, = 7'e*-*’, and integrating J, by parts§ 


we have 2-1/2 
 — |2p2—2vz Wes. ea 1 \ 
f(2) = we 23I2, [14+ o(;)| 


uniformly in @. Hence at a zero 
e2 — 7- [%g2v2-1/29-3/22-16 ] 1 O(r-1)}, 
Hence by Hardy’s method 
2a = 2r¥2 cos 36—} log(8ze)—log r + O(r-), 
2y = are sin 30—0+ O(r-). 
Putting r= y(1+e), 6 = 4m(1—e), we have (9.11). 

9.2. It may be observed that, although a function of the form 
(1.11) may have a few real zeros, Titchmarsh’s|| theorems limit the 
real parts of the main stream of zeros very strictly. For > |cos6@, |/r,, 
is convergent. This prohibits x, = n/logn except for a comparatively 
small proportion of zeros, although n(r) ~ 2r/7 permits 

Y, = +nr+n/logn for all n > 2. 
This limitation of the real parts of the zeros is a characteristic of 
functions which are the sum of exponentials, and by definition f(z) 
is a limiting case of such a function. 

+ Z, Theorem XIII. 

t E. Goursat, Cours d’ Analyse, 3rd ed., vol. i, p. 295, Ex. 23. 

§ Compare G. H. Hardy, Proc. London Math. Soc. (2 (1904), 405. 


2) 2 
E. C. Titchmarsh, Proc. London Math. Soc. (2) 25 (1926, Theorems I 
and IV. 








SOME APPLICATIONS OF GENERATING FUNCTIONS 
TO NORMAL FREQUENCY 
By A. C. AITKEN 
[Received 20 February 1931] 
1. Generating functions of frequency and moments 

In many problems of mathematical statistics the moment-generating 
function is easier to handle than either the frequency function or the 
frequency-generating function. Indeed, the form of the frequency 
function can often be deduced by inverting the moment-generating 
function. The conditions for the validity of this inversion, and for 
the convergence of moment-generating functions, have been very 
fully discussed in an important paper* by T. Kameda. No com- 
plications arise in the cases treated below, the frequency function 
being the exponential of a negative definite quadratic form, the 
character of which may be preserved by a suitable choice of para- 

meters « in forming the generating function. 
The intention of the present paper is purely expository, to give 
some examples of the solution by moment-generating functions of 
some important problems of normal frequency. No new results are 


derived, but the systematic use of a determinantal lemma, the 


extension of which is interesting on its own account, seems to have 
certain advantages. 

The generating functions are defined} thus. If 4(2,, a9,...,x,), or 
briefly ¢, is the frequency function of » variables, then 


> >d Dd...b. Gp... 2, or [ [ [ wo D.C? 42" daday...dx,, 


according as the variables are discrete or continuous, is the generating 
function of frequency and, with a substitution ¢ = e%, 
>> >d---¢-exp(> ax), or || | ...d-exp(> ax) dada, ...dzx,, 
is the moment-generating function, since a typical term in «, 
ojtadz... ain/(i,!t9! ...%,!), 
has for coefficient 
= sty mle at 
N46, Agen tg z z > oe n®, 
or the analogous integral. 


* T. Kameda, Journ. Fac. Sci. Imp. Univ. Tokyo, (1) 1 (1925), pp. 1-62, (38). 
+ H. E. Soper, Frequency Arrays, Cambridge, 1922, pp. 8-13. 
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By the multiplication theorem of frequency, the moment-generat- 
ing function of the sum of several independent variables, or of 
independent values of different simultaneous variables, is the con- 
tinued product of the several moment-generating functions. 

To change the origin of measurement of x, from zero to h, implies 
the writing of ¢7-~’ for t7, in each term of the frequency-generating 
function; hence if F(«) is a moment-generating function, the corre- 
sponding moment-generating function about the mean (%,, Z, 


is exp(—} a&). F(a). 


2. Quadratic forms and a special determinant 
From an abstract standpoint the theory of normal distribution in 
several correlated variables is the theory of the linear transformation 
of a positive definite quadratic form, such as 


= C(x, 2), 


in so far as the properties of exp{—}C(x,x)}, its integral and its 
moments are affected. Certain standard results are therefore re- 
quired. First, 


4) 


, .. exp{—4C(z, x)} dx,da,...dx,, = (2m)"C-4, (2.1) 


— 
where C = |c,| is the determinant of the form C(x, x). 
This well-known result follows at once from the fact that when 
v,«) is transformed into a sum of squares > £* by a transformation 


: , 2 » ae acobis on it 
> Aj; &, then |h;; c;;|-1; hence the Jacobian |h,;| = |c,;|-*. 
J 


Next, pete: v,x+v) = C(: +2 > ax+C-\v, v), (2.2) 
where C-}( =} claja;, and c= C;,/C. Thus C-\(a, «) is the 
reciprocal pads i atic form, the a’s being defined below. 
C(a+v,a+v)=> C;;(%;+¥;)(4;+¥,) 
— > Cj,X Xj; +2 > x; > Cj + > CijViV; 
= O(x,x)+2 >} ax+C-"(a, a), 
where «; = > ¢;,v;, for j = 1, 2,..., , and so v4; = > ca;. 
The point of this is that just as in elementary algebra we ‘complete 


For 


the square’ of ca®+2aa% by adding c-'a?, so a quadratic form plus 
linear terms, C(x, x) +2 > ax, may be ‘completed’ by adding C~\(a, a). 
Lastly, the values of two related types of determinant will be 
required, 
K2 





132 A. C. AITKEN 
a+b b B aR A+B B B 


| 


b atbdb b...| 23 1 B A+B B 
(2.3) anc B B At 


1 


b b a+tb ata 

; ha 

In (2.3) the diagonal element is always a+-b, and every other ele- 
ment is b. The operations ‘row 1 —row 2, row 2 — row 3....,’ and 
then ‘col n+ col 1+ col 2+...+ col n—1’ reduce it to a form in 
which its value* is seen to be a”-1(a-+-nb). In (2.4) the A and B repre- 
sent square m-by-m blocksof elements, the block A + B being obtained 
by adding corresponding elements of A and B. The procedure for 
evaluation is exactly similar to that for (2.3), rows and columns of 
blocks being subtracted or added like elements. Laplacian expansion 
then gives the value of (2.4), if of order mn, as |A|"~1. |(A+nB)|. 


3. Normal frequency in many variables 

We proceed to investigate the simultaneous frequency-function of 
variables 2, X,...,Z,, where each x;, measured as a deviation from 
its mean value, arises as the sum of NV elements 5z;, each simultaneous 
set of elements being subject to a law of frequency. The law is taken 
to be of such a kind that the mean absolute value of |5x;| is of order 
N-1, the mean values of (5x,)? and 6x,5x; of order N-*, while the cubic 
and later elementary moments are of higher order. These conditions 
are satisfied by a very wide class of laws of frequency, but not by 
all; for example, in the Poisson case, where the probability of non- 
zero deviations is itself of order N~-!, the elementary moments are 
of higher infinitesimal order than the above. 

The moment-generating function of x,, %,..., 2, will be the con- 
tinued product of N elementary moment-generating functions, and so 
the second or quadratic moments 7;;0,0; are of order N-1, the cubic 
and later moments of higher order. The moment-generating function 


| e : ip! eee , as 
l — > Vj FT jXjOX; 2!-+ ioe (15; Ss i. Yj o- riz). 
is therefore replaced, to the same order in NV, by exp{}R(ca, oa)}, 
where R = |r;;|; and the continuous frequency-function ¢ which has 


ij 


this moment-generating function is given by 


exp{} R(ox, ox)} = [ [ [ ..6.exp(> ax) dx dxq...dx,. 


—2 


* The determinant is well known and is usually evaluated indirectly by con- 
siderations which, however, do not apply at once to (2.4). 
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The solution of this equation, by (2.1) and (2.2), is 


$ = (27)-!"| R| 0,0, ....0, exp{—4}R-(2x/0, x/a)}, 


and this is the Edgeworth-Pearson function of normal correlation. 
It arises also if the 2,, %,...,%, are linearly correlated* through 
being each of them linear functions of a larger number of uncor- 
related variables €,, &,..., &, normally distributed. 


8 


4. Distribution of normal variance 


For the next example we may investigate the distribution of s?, 
the variance or squared standard deviation of x, as calculated from 
samples of N values drawn from an infinite normal stock with 
frequency-function (27)-te-#”". Following Wishart? we shall denote 
s° by a). 

A statistical parameter @ will in general itself possess} a moment- 
generating function 


F(0, x) = [ $(4)e% dd, 


where, if 6 is calculated from independent sample values 2, %9,..., ,,, 
the function 4(@) will represent the compound frequency of the 
sample values, and the integral will be N-ple. Thus in the present 
case the variance about the mean of sample is 

a,, = > 2?/N—(> z)?/N?, 
and so F(a,,, a) 


= (2Qer)-14 [ [ [ ...exp(—} > a*)exp{a > 2?/N—a(> x)?/N*}da,dxq...dxy. 


The determinant of the quadratic form under the exponential is 


of the type (2.3), with 1—2a/N for a, 2a/N? for b. Its value is thus 
found at once to be (l1—2a/N)*-1, and so, by (2.1), we have 
F(a,,, «) = (1—20/N)40-9 = 14(N—1)o/N +... 
The coefficient of « puts in evidence the familiar result that the 
mean of s? for samples of N is (N—1)o?/N, where o? is the variance 
of the sampled population. Also the moment-generating function of 


* M. J. van Uven, Proc. Amsterdam Acad. 16 (1214) pp. 8--13. 

+ J. Wishart, Biometrika, XXa (1928) pp. 32-52. Also V. Romanovsky, 
Comptes Rendus, 180 (1925), 1897; Metron, 5 (1925), 27, 32, 35-6. 

t Kameda, loc. cit. 
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s* about its own mean is then 
(1—2a/N)-#4-) . exp{—(N—1)a/N}, 
from which its own variance and higher moments are easily cal- 
culated. 


5. Product-moment in bivariate normal sample 
We next consider N sample drawings from an infinite normal 
bivariate population with frequency-function 
1 
(27r)-1(1—p?) texp{ — 1 2 , («?—2pary+y?)}, 
oe 


and suppose the product-moment 78,85, Or @,., to be calculated about 
the means as 4 72 

Ayo = > ry/N—(> x)(> y)/N?. 
Thus the moment-generating function, say F(a,., «), is the 2N-ple 
integral of 


1 
> 


i ae a oe eal le cal \ 42 
(277)-¥(1—p?)-!Nexp a> ¢ 2pry+y")) x 


| 
xexp{a > ay/N—a(> x)(> y)/N?}. 
The determinant of the quadratic form in the integrand now 
belongs to the type (2.4), where the arrays A, B are seen to be 


0 a/N? 


Hence 


—}(N—-1) 


F (ay, «) ] —p—7(1—p") 


xX 
Sets | 1 
p v! p”) 


= {1—2pa N —a?(1—p?) N?} (N-1), 

The mean value of rs,s, as calculated from samples of N is seen to 
be (N—1)po, o,/N, and the moment-generating function about this 
mean is 

exp{—(N—1)pa/N}. {1—2pa/N —a?(1—p?)/N?}-#4-9, 


6. Simultaneous quadratic moments in sample 


The two preceding examples indicate a general process and result. 
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Samples of NV are taken from a multivariate normal population with 
frequency-function 

(277)*" |p;,;|-texp{—3 R(x, x)}, 
n 
where R(x,x)= ¥ p,x;x,, and R-\(x,x) is the reciprocal form. 
i,j=1 


Quadratic moments r;; 8; 8;, or a;;, are calculated by means of 


a;;= >’ x,x;/N—(>’ x,)( >’ x;)/N?, 

where >’ is used to denote summation over the N sample values of 
particular variables, and >} is reserved for summation over the n 
different variables. 

The simultaneous moment-generating function of sample quadratic 
moments, F(a, «), is thus the nN-ple integral of 
(2n)-" |p,,|-Nexp{—} 5’ R(x, x)} x 

xexp{N— >” ¥ «2%; 2;—N-* > a,,(>’ x;)(d’ 2;)}- 

Once again the determinant of the quadratic form involved, of 
order nN, is of type (2.4), the arrays A, B being 
Piles Gye s ss i 2041 %2 


B 9 


ye] %2 “%2- + «I 


oe 1 
A =" i= yo 2g. - «|? 


If the array of «;;’s be denoted by [«;,], with the convention that 
diagonal elements are doubled, then by (2.4) 

F(a, x) = |pjj|-* |[p? —N-*a;]| #4 —-? |p |-4 
= |(1—N+ [pis las])| sae, 

The above denotes a power of a determinant in which the elements 
are those of the product —[p;;|[«;;], multiplied by N~ and increased 
by unity in the diagonal. The result, allowance being made for 
difference in notation, is that first given* in 1928 by Dr. J. Wishart, 
who obtained it as a deduction from his exact form for the simul- 
taneous frequency distribution of these quadratic moments. 

The corresponding distributions, not about the mean of sample 
but about the mean of the population, are obtained at once by sup- 
pressing the terms in 1/N? in the formulae for calculating the 
moments. The effect is to replace the elements b or arrays B in the 
determinants by zeros, so that the power of the resulting determinant 
appears finally as —}N instead of —}(N—1). 


* Wishart, loc. cit. 








SIMULTANEOUS QUADRATIC AND LINEAR 
REPRESENTATION 
By GORDON PALL 
[Received 20 January 1931] 
WE investigate the solvability in integers 2,, ..., x, of the simultaneous 
equations — y O50, 2; 7 | ae (1) 
b= > b,x; 
where f and / are respectively quadratic and linear forms with integral 
coefficients, and 
Aj; = Ay; D= |a;;| #9. 
The special equations 
b=a,+...p2 


- 2 _| . 
a=AUMA7---TIXs; 


9 
8 
) 


are studied in detail (sections 2, 4, and 5). 

It is seen that there exists a (1,1) correspondence between the 
integral solutions (x, 
of an equation k(Ka—Db®) = (yo,---, ¥,) (4) 
and a system S of linear congruences in the y, and b. Here ¢ is 
a quadratic form with integral coefficients, k is a constant, and 

K = F(b,,..., 5), (5) 
where F is the adjoint of f, is supposed not zero. 

The system S can be reduced to simplest terms, and in special 
cases can be entirely eliminated by using the signs of the y, or other 
transformations of . 

In section 3 we shall find a method of obtaining a set of sufficient 
conditions for the solvability of (1) in integers x; >0, when f is 
positive definite. 

An interesting way of expressing the above-mentioned (1, 1) corre- 
spondence is illustrated in Theorem 2. 

In section 5 we shall, for any s > 8, determine all pairs (a, b) of 
integers such that equations (3) are solvable in integers 2;,. 

1. Let F be the adjoint form of f. Then 

det(ra,;—b,b;) = 7°D—7F*F(b, 
Hence the determinant of the form 


$lat,,...,%,) = Kf—DP 
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is zero. By rational operations we can express it in the form 
kb = 3 ¢,X3=WX,,..., Xy) (8) 
Xp = ChyXytCpgVot--- +p, (h = 2....,8), 
where k and all the c’s are integers, k > 0. 
By (7), (8), and (1,), 
kKf = kD(> b,x;)?+0(D ¢o,x;)?+...+e,(¥ 2) 
The Hessian being an invariant of weight two, 
(kKK)8D = kDe,¢3...c,A”, 
where ¢,;=b,; (j= 1,...,8), A= |e,|. If KAO, A and ¢,....,¢, 
not zero. Hence the equations 
b,x,+...+ba, =b, 
Cyj% +... +6, %,=Y, (h=2,...,8) 
have the unique solutions 
= (Cyjb+ Co; Yet... +C,;y,)/A, 
where C;; is, naturally, the cofactor of ¢;;. 
THEOREM 1. Let f and | be as in the introductory paragraph, and let 
b,) 40. We can set up identities (8) by simple operations. 
For any such identity there is a (1,1) correspondence between the sets 
of integers x; satisfying (1) and the sets of integers y,, satisfying (4) 
together with a system S of linear congruences involving the y, and b. 
Evidently S is any system equivalent, under the hypothesis that 
the y;, satisfy (4), to 


C1 jb+ Co; Yet... +C,;¥, = 0 (mod A) 


2. For (3), one expression (8) is 
(s—1)!(sa—b?) => ia (h - i) Xj, (14) 
(} — 


where Xp = (h—1) 25h —Vp—nsq— +» —Xg- 

As generally for (1), the identity (14) for (3), obtained by com- 
pleting squares and preserving the integrality of coefficients by 
adequate multipliers, can be greatly simplified by using relations 
between the X,. For example, the identity for s = 4, namely, 

3(4a—b?) = X34+-2X24 6X3, 
where X,= X,+3Y, Y =2,—2z,, readily reduces to 
4a—b? = (x, +2,—2%3—2,)?+ 2(a%3—2,)?+ 2(4,—29)”, 
and hence, by way of the almost trivial identity 


262-4 2n? = (E-+n)*-+(E—a) 
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to the classical identity 
4a—b? = (x, +2,—%3—2,)*?+ (4, —%.+%3—2,)?+(x 
In a similar way, for s = 3, 5, 6 we get 
4(3a—b?) = (7, —2,)?+ (2, —2_)(%.—2%3)+(%_2—2)*; 


4(5a—b?) = (47, —x7,.—X3—2%,—2;)*+ 5(2,+2%3—2,—25)° 


2(6a—b?) = (22,+ 24,—xX3—%,—2,—2%4)?+3(7g+2,—2,—2Xe)?+ 
+ 6(x,—2_)?+ 6(x,—2,)*+ 6(x;—2)?. 


A repeated use of (16) gives a simple identity for 





s(v?+...+27) (e == 3) 
as a sum of 2” squares (x,+¢,%.++...+-e,x,)?, where the e, are +1. 
Let N,(a,b)-denote the number of representations of the integer 
a as a sum of s squares the sum of whose bases is 6; that is, the 
number of sets (2,,...,%,) satisfying (3). 
THEOREM 2. Let a and b be integers. Then 
(A) N,(a, b) = w,N(6a— 2b? = a?-+-3y?) 
= w,N (3a—b? = 227+ 2xy+ 2y?); 
(B) N, (a,b) = w,N (4a—b? = 2? +-y?+-2?) 
= w,N (4a—b? = x?-+ 2y?-+ 227); 
(C) N,(a, b) = ws{N(5a—b? = a? + 5y?+ 527+ 5w?)+ 


+N (20a—4b? = 2*+ 5y?+ 527+ 5w*)} 


where w,=we=1 (3)b), W,=—w,=1 (a=b=0), 
=4 (3/6); =} (a=b=1), 
w, = 4 (5|b,a=b), =0 (afb); 


- 4 (5/b,a =b), 
= 0 (a #b), mod 2. 

A. Let s=3. Equations (3) are algebraically equivalent to (18) 
and (3). Hence there is a (1,1) correspondence between the sets of 
integers x, satisfying (3) and the sets of integers y, satisfying 

}(3a—b*) = ¥3+Y2Y¥stY3, 46 =Y2—Ys (mod 3). (21) 
Now by (21,), 6? = (y.—ys3)" (mod 3). Hence if 3|b, (21,) holds for 
every set of y,, of (21,); if 3/6 it holds for precisely half of these 
sets since but one of y,—Y3, Y3—Y2 is congruent to b (mod3). Hence 
we have the value of wy. But w, = we, since 
N(4n = x?+3y?) = N(n = 2?+2ry+y?"). 
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B. Let s=4. On using (17) in place of (18), we see that there is 

a (1,1) correspondence between the x; of (3) and the y,, of 
4a—b? = y3+y3+ 93, YotYyzt+y,+b = 0 (mod 4). (22) 

For later use we note that the correspondence is given by 
® = HO+Y2+Yst+Ya), t= }(6+Y2—Y3—Ya), ete. (23) 
If a#b (mod 2), (3) are of course not solvable. If a,b are odd, 
(22,) implies 2, y3, y, odd, whence (22,) holds for precisely half the 
solutions of (22,), just one of b+y,+y,-+y, being divisible by 4. 
Now c? = 2c (mod 8) if ¢ is even. Hence, if a, b are even, (22,) implies 


—2b = 2y.+ 2y,+2y, (mod 8) 


and hence (22,). Thus we have ws. But w, = ws, since 
2 


N(n = x? + 2y?+ 22?) = N(n = 2? + y?+2?) if n = 0,3, 4, (mod 8). 
C. Let s=5. Using (19), we find a (1,1) correspondence between 
the solutions x; of (3) and y, of 
4(5a —b?) = y3+-5y3 + 5yj + 5y3, (24) 
Yo = 6b (mod 5), yot+y3+y,+y; = 0 (mod 4). 


The rest of the argument is left to the réader. 


3. Inequalities; solution in integers x; > 0 
Let f and lJ be as in section 1 and let f be positive definite. Then 
K >0, D> 0, ¢g,...,¢, > 0, and so 
Kf > Di? for all real 2;. (25) 
As a corollary, if f is positive the inequality Ka > Db? is necessary 
for the solvability of (1) in real 2;. Suppose 6; >0 (¢= 1, 
f=a,7}+...+a,2?. (26) 
The assumption b > 0, x; << —1, for some 7, implies 
(b6+6;)? <7,(a—a,), 7; = F(b,,...,b;_,, 0, 5;,1,---, 6,)/D, 
by use of (25). Hence, if f and / are of the above form, the x; of 
every solution of (1) in integers are all non-negative if 
b>0, ‘ , 8). (27) 
To get a set of sufficient conditions for solvability of (1) in integers 
x; 20, merely adjoin (27) to a set of necessary and sufficient con- 
ditions for solvability im integers 2;. That (27) may sometimes be 
improved on is illustrated in the next section. 
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4. Conditions for solvability of (3) in non-negative integers 
By Theorem 2, (A), equations 
a = 27+ -y?+-2?, b=a+y+z 
are solvable in integers 2, y, z, if and only if 
}(3a—b?) is represented in €*+-37?. 
By section 3, each of x, y, z >0, if b > 0, b?-+2b+-3 > 2a. 
THEOREM 3. If s = 4 the conditions 
a = b (mod 2), 4a—b* =a sum of three squares, 
are necessary and sufficient for the solvability of (3) in integers x; No 
x; of any solution is negative if 
b>0, 624+2b4+4>32. (31) 
A solution in integers x; > 0 exists if (30), (31,), and 
a,b are even, 362+ 8b+ 16 > 8a. (32) 

Only the last statement is new.* Its proof is clear from (22) and 
(23). When a, b are even, the signs of the y, are all at our disposal, 
and, given a representation (22,), we can suppose Y, < Y3 < —¥Y4, 
y, 0. Then the x; of (23) are all greater than —1 if 

—Y2—Ys—Y4 < 5+4, 
and hence, since 
(—Y2—Ys—Ya)® < (|Yo|+ lys|)? < 2y3+y3+Yi) = 2(4a—6?), 
if 2(4a—b?) < (b+-4)?, which is (32,). 

THEOREM 4. Jf s=5, 6, or 7, equations (3) are solvable in integers 
%, of and only if a =b (mod 2), sa > b?. (33) 
They are solvable in integers x; > 0, if also 

b>0, b? > 3a—5. (34) 


I. Proof for s=5. We use the fact that an integer is a sum of 
three squares, if and only if it is positive or zero and not of the form 
4"(8v+7), to show that 5a—b? is represented in 2?+-5(y?+-2?+w?), 
if and only if 5a >b?. Hence the first part of Theorem 4 for s = 5 


follows from Theorem 2, (C). 


* Cf. Cauchy, @uvres, (2), vol. 6, pp. 320-53; Legendre, Théorie des 
Nombres, ed. 3, vol. ii, Nos. 624-30. Many elaborations and extensions of 
Theorem 3, beyond those of the present section, are to be found in the writer’s 
paper ‘On sums of two or four values of a quadratic function’, to appear 


shortly in some American journal. 
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If (6+-1)? > 4(a—1), b > 0, then all x; are positive or zero by the 
preceding section. Hence assume 

3a—5 <b?<4(a—1), b>0. (35) 
If a, b are odd, take x; = 0; if a, b are even, take x; = 1. By (35) all 
conditions of Theorem 3 are satisfied with a—x?, b—x; in place of a, b. 

II. Given that equations (3) are solvable in integers when (33) 
holds, [non-negative integers, if (34) holds too] for s= S—1 >5, 
we wish to examine its like solvability for s = S, if 

(S—l)ja <b? < Sa, (b> 0). (36) 
We have to seek a non-negative integer 2 such that 
d(x) > 0, d(x) = (S—1)(a—x*?)—(b—2z)?. (37) 
This must hold for x the integer nearest b/S, which is the turning- 
point of the parabola z = ¢(y). The only case of doubt is therefore 
¢(+4+6/S) <0, whence a and b satisfy 
Sa > b? > Sa—}8S?/(S—1). (38) 

If S =6 or 7, (38) becomes Sa = b?, since Sa—1.is not a square. 
Then b = Sx, a= Sx, and (3) are solvable with z;=x. This com- 
pletes the proof of Theorem 4. 

If S = 8, (38) reduces to b®>= 8a. Then b=0(mod4). If b= 82, 
then a = 8x? and (3) are solvable as before. But if b= 4y (y odd), 
then a = 2y?, and equations (3) for s = 8 are not solvable in integers 
«;. For (37,) holds for no integer 2 in view of ¢(4y+ 3) = —2. 

5. Case s> 8; all (a, b) satisfying (3) 

Let (3,) denote equations (3) for s = 8. 


LemMA 1. Let s>8. Let dg, by be integers of the same parity for 
which equations (3,) are (i) solvable, (ii) not solvable, in integers 2;. 


W rite r = 8a,—bj, p = least absolute residue of by (mods). (39) 


Then equations (3,) are (i) solvable, (ii) not solvable, in integers x; for 
every pair (a, b) of integers satisfying (33,) and 
sa—b* =r, b =p (mods). (40) 
For any integers r, p, 8, (s >8), let (r, p), denote the set of all 
pairs (a, b) of integers of the same parity satisfying (40). This will 
usually be a null class (compare Lemma 2 following). The set of all 
(a+2*, b+2), for all (a, b) of (r, p), and all integers x, coincides with 
the set of all (a’, b’) of like parity satisfying 
(s+1)a’—b" =7’, b’ = p’ (mods-+1) (41) 
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’ for any integer v, where 
r’ =r+{r+(sv—p)}/s, — p’ = p+. 
To see this we calculate (s+ 1)(a+a?)—(b+<a)? with x = v+(b—p)/s. 
The sets (r’,p’),,, are said to derive from (r, p),. 

Now all the pairs (a,b) of integers of like parity can be segregated 
into sets of the form (r,p),. Call a set (r, p), solvable if (3,) are solvable 
in integers 2; for all the elements (a,b) of (r,p),; unsolvable if (3,) are 
unsolvable for all (a,b) of (r,p),. For s = 8 every set (r,p), is either 
solvable or unsolvable; indeed (0,0), and all (r,p), with r >0 are 
solvable, and (0, +4), and all (r,p), with r <0 are unsolvable. 

Every non-null set (7’,p’),,, derives from an infinity of sets (r,p),, 
all obtained by solving (42) for r and p. Every derived set of a 
solvable set is solvable. From the nature of the transformation 

(a’ = a+2?, b’ = b+2; all x) 
it is now clear that, if every set (7, p), is either solvable or unsolvable, 
the same is true of every set (r’,p’),,,. From the fact for s = 8, every 
set (r,p),, 8 > 8, is solvable or unsolvable. This is the lemma. 

Lemna 2. A set (r,p), 1s non-null, of and only if 

r+ p?—ps = 0 (mod 2s). (43) 

For this is a necessary and sufficient condition for the existence 

of an integer ¢ such that 


9 


(p-+ts) 


8 


- : , ‘ 
is an integer of the same parity as p+ts. 


Lemma 3. Let r, p, 8 be integers satisfying (43), s > 8, and |p| < $s. 
Then the following three propositions are equivalent: 
the set (r,p), contains an element (a,b) with a <b; 
the set (r,p), 1s unsolvable; 
r< |p|s—p?. 
It is trivial that (44) implies (45). For (3,) imply 
a> > |z,|>6. 
It is trivial that (46) implies (44). For let b= p if p>0, b=p+s 
if p <0. Then a=c or c+2p+s, and a <b, where 
c = (r+p?)/s. (47) 


We prove by induction on s that (45) implies (46). We know this 
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to be true if s=8. Assume for a given s >8 that the unsolvable 
non-null sets are precisely those satisfying 

(43), (46), and |p| < 4s. (48) 
The sets (r,p), from which (r’, p’),,, derives are, as we see by inverting 


etl 
9 y] yr , 
(42), given by lee r’ +4(s-+1)o—p'}? 
s+1 
where v = 0, +1, +2,.... The set (r’,p’),,, is unsolvable, if and only 
if (46) holds with r in (49) and with +, in place of p, for every v. 
Here 7,, denotes the least absolute residue of p’—v modulo s. Hence 
(r’, p’),,, unsolvable implies 
r’ —{(s+1)v—p’}?/s < (s+-1)|7,,| —(8+1)72/s (50,) 
for every v. If |p’| < 4s, (50p) is 
r’ <(8s+1)/p’|—p”. (51) 
If p’ = }(s+-1), (50,) is equivalent to (51). The induction is complete. 


> 


(49) 


As a consequence of these lemmas we have 
THEOREM 5. Lets >8. For each of p= 3, 4,..., [48] determine all 
integers r > 0 such that 


™~ 


rp? = r+ 
. tas 


p (mod 2). (52) 
8 


For any such p, r and any a, b satisfying 
b = +p (mods), sa—b? =r, (53) 
the equations (3) are not solvable in integers x;. They are so solvable 
for every other pair a, b satisfying (33). 
The largest r is [}(s?—8s)]. All values (p,7) for some small values 


s follow: 














INTEGRALS EXPRESSING PRODUCTS OF 
BESSEL’S FUNCTIONS 

By T. W. CHAUNDY 

[Received 20 March 1931] 

1. Mellin’s principle 
THERE is a principle due to Mellin* which enables us to obtain many 
of the standard integrals in the theory of Bessel’s functions (and 
elsewhere) in a systematic and comparatively simple fashion. Since 
I cannot find reference to it in Watson’s Bessel’s Functions, it may be 


that this method deserves to be better known. 
The integrals ome by Mellin are essentially of the type 


-{ (Z)yo%, (1) 


taken either between Limite ecole of z or round a closed 





y= 


contour. 
I write 5 for the operator zd/dz. Then 


by = ‘| p’ (Fu 2B = dt 
. —|4(F)vo] | 4(7) {tp’(t)} = (2) 


on integration by parts. It is convenient to write { for z/t, to denote 
the operators in (,¢ corresponding to 6 by 8,, 8,, and to regard these 
operators as essentially partial differential operators operating only 
on the functions of f,¢ respectively to which they are attached. We 


may therefore write (2) as 
. . y 
5 | Hosen F = | HORM) F [HHO 


where the brackets [ | indicate terms dais at the limits. Since 
54(C) is 8; $(2), we have by repeated applications of the above process 
~ r dt 1) n n—1S | n 

5» f omy = | $CpBr oe) F— [0p 8p8, +... EY BOMO)) 


and Saiilien more generally, if f( ) denote any polynomial with 


“3 


constant coefficients 


8) { some S = [ sopennn SF — [ML sanun]. 


8:—3, 


* Acta Soc. Fenn. 21 (1896), chiefly § 9 








ee 
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Tf f,, fos 1s Ja ave any four such polynomials, we can then write 


f,(8)fs(8)—zg,(8)go(8)} | 4(0) 
a lg Ri i ~ = 


—| {neo peP- = — gmmgy 5) O80) gow]. (4) 
ee ee 


The argument presupposes thet the differentiation across the integral- 
sign is permissible. In the particular applications which follow this 
can generally be seen to be so without special discussion. 

If now ¢, # satisfy the respective differential equations 


(f:(8)—2"91(3)}P(z) = 9, {fa(3)—2g2(8)}¥(2)= 0, (5a) 
the integral on the right of (4) disappears, and we have only to secure 
the evanescence of the ‘terms at the limits’ in order to secure that 


y= | dict) 


be a solution of the differential equation 
8) f2(8)—2"91(8)g2(5)}y = 9. (5B) 

If the integral is a contour-integral, the terms at the limits vanish 
automatically on completion of the contour, provided that the expres- 
sion in square brackets is one-valued. Our task is then to see that 
the contour-integral is neither infinite nor zero. If on the other hand 
the integral is a rectilinear integral, we must choose limits of integra- 
tion such that the expression in brackets takes the same value 
(invariably zero, as it happens) at each limit. Here I shall consider 
only the rectilinear integral. 

In the terms at the limits we can write 


fb POLED goog 
= 


aie 0 mg (By )b(CWh(t 


—- {@ 5,22 (3¢-+m) —fo(9,) 
= ("g,(3;) 8;-+m—8, (f(t). 


Thus the expression which must vanish at the limits may be written as 
cm 5,) (26 ¢+m) —fo( (8,) —_ fm J2 (8z)— J2(5,) | t). 6 
Saal 0) 5¢:-+m—8, 5;—5, jee We) (6) 
L 
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We may accordingly enunciate Mellin’s principle as follows: 
The differential equation 
fi(5) fo(8)y = 2"91(8)g2(9)y 
is satisfied by the definite integral 
y= [ ooo S 


. 


if ¢, % are solutions of the respective equations 


f,(5)¢ = 29, (8)4, Sa(3)yb = 29 2(3)x, 


and if the limits are chosen to make the expression (6) vanish. 


2. Prolegomena 
In illustrating this principle by its applications to Bessel’s func- 
tions, I confine myself to those integrals which express products of 
pairs of Bessel’s functions. Except in limiting cases I suppose these 
functions to be of non-integral orders yp, v. 
The differential equations satisfied by products of pairs of Bessel’s 
functions of the same kind are orem 


{82—(u+-v)}{32@—(w—v) 9} ny Sn = —42%(8+1)(8+2)Tay Tay (7) 
{3*—(-+) oe p— Lay Loy =422°E+WE-+2)L ny lay. (8) 
If p=», rie se bsit: become fied 
82(52—4v2)y = —422(8+1)(8-+2)y, (7A) 
82(52—4v2)y = 422(5+-1)(8-+-2)y. (8a) 
Evidently (7A) has the solutions 
< . in 
and the ‘missing’ fourth solution is supplied by 
£. tO Ody 
Ov Ov 


If we remove this fourth solution, we get the equation 
6(82—4y?)y = —42z7(5+-1)y, (7B) 
which has just the three solutions 
J?, : d*.. 
By change of sign of 2? we obtain similar results for (84) to which 
corresponds the ‘reduced’ equation 


8(82?—4y2)y = 42°(8+ 1)y. (8B) 


* Watson, loc. cit., p. 146, 5.4 (3). 
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The chief analysis of equations (7), (8) turns on the pair of dif- 
ferential equations 
(8°—p?)F., = —2%(5+1)(8+2)E.,, (9) 
(8°@—p2)F.,, = 2°(8+1)(8+2)F.,, (10) 
of which the solutions are given by 


1 2)-+1)4 
E,(2) = vag EN (11) 
THE V(1—2?)+1\# 
v(1—z?) | z i 
At the origin, we have, if » is not an integer, 
E (2), F(z) = (32) *+O(2 *), 
and at infinity E,,(z) = 21+ O(z-). 


I add for reference that 


Fi, 


(12) 





F(z) = 


E_(cosech uw) = tanh u e#“ 
p 2 


F,,(sech u) = cosh u er, F,,(sec u) = —icotu et, 


3. Nicholson’s integral for K,K, 
We first analyse (8) into the equations 
{8° (uv) (2) = 424 (2) 
(5*—( Aas: pryaieogiyor 
Taking ¢= K,,,,(2z), $= , we get the Mellin integral 


2) 
jn) 0 
The terms at the limits are now 
407(1—12)(8;+8,+:3)—D}K,,,,(20)F, 
These vanish as t>0-+-, if R(z) > 0. Writing t = sech uw, we therefore 
have as a solution of (8) 


i 2) 
( K ,.,(2z cosh uje#—” du, 
2) 


or more conveniently, 


[ K,.,(2z cosh u)cosh(u—v)u du. 


0 


pty 


L2 
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Now cosh wu > 1 and the factor K., ,, therefore assures at R(z) = +00 
a zero of order less than that of any 
e-(2-O)R2) (e +> 0). 
The general solution of (8) is a sum of constant multiples of the 
products , , - 
K,K,, £.i,. ee Egle 
We obtain a zero of exponential order, only if we take the first of 
the products in isolation. We therefore write 
Ms 
z)K(z)=C | K,,.,(2z cosh u)cosh(u—v)u du. (17) 
0 
The value of C is found to be 2 by identifying, on the two sides of 
(17), the dominant terms at z= oo. We thus reach the final formula: 


ul 


K, (z)K,(z) = 2 | K ,4(22 cosh u)cosh(u—v)u du, (18) 
0 
where R(z) > 0 and p, v are unrestricted. This is Nicholson’s formula, 
Watson (loc. cit), p. 440, 13.72 (1). 
4. Integrals for /,,/, and J,,/, 

Still contemplating the analysis of equation (8) into equations 
(15), (16), we can utilize J,,,,(2z), the other solution of (15). This 
leads to the Mellin integral 


| a a 
“sil | Iyo(F) F-y(t) 7? Say, 


= | I,4,(2zc0s u)e—™ du, 


if we write t= secu. The terms at the limits now become 
2z7{sin 2ud;-++ 4 cot u(1+-sin*u)—i(u—v)cos*uje™ "I, . (2), 
where € is 2zcosu. These vanish at u= +47 if R(ut+v+1)>0. 
With this condition we have the new solution of (8): 
hor 
I ,+,(2z cos ue“ du, 
x os 
or more conveniently 
br 
| Z,,+»(22 cos w)cos(u—v)u du. 
0 


We can expand this integral in powers of z, the indices which occur 
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being the set 4+v-+2r, where r is a positive integer or zero. The 
general solution of (8) is a sum of constant multiples of 


—, i 2 I_,1 


pow —p*y? —p* -9* 
If 2u, 2v are not integers, we obtain a sum of powers z#+”+*" only, 
only from the first of the products in isolation. We therefore write 
ia 
LI,=C i] I 


0 


ptv(2z Cos w)cos(u—v)u du, (19) 
extending the results to integer values of 2u, 2v by continuity. C is 
found to be 2/7 by comparing, in (19), the leading terms at z= 0. 
We thus get at length 
bar 
srl I, = | I 
0 


p+v(22 COs u)cos(u—v)u du, 
if R(utv+1)>0. 
If we write » = —v, we deduce 
40 
inl I_,= ( I,(2z cos u)cos 2vu du. 

0 
Again, from (20) we have 
ta 


elt @ See eee fs" 3 en 
. | Beye "Op “) J &u+r) p+r(2z COS U)cos(u—v)u 
0 


If again we put » = —v, we deduce 


47 
0 é . 
in(,- 1-121, = | K,(2zcosu)cos 2vu du. (22) 
. ov ov 


We may similarly analyse equation (7) into 

{8°?—(u+v)*}6(z) = —424(z) 
and equation (16). If we take the solution J,,,,(2z) of (23) and pro- 
ceed as we have just done with the sblation I ,+y(22) of (15), we 
arrive at the formula, —e to (20), 


bad, J, 2 = Ly 
0 
where R(u+v+1)>0. From this we may similarly deduce 
tn 
bad, J_, = | J,(2z cos u)cos 2vu du 
0 


p+v(22 COS u)cos(u—v)u du, (24) 
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tar 
and +(J., I-Ie I) = | Y,(2z cos u)cos 2vu du. (26) 
0 


Of the foregoing formulae, (20) and (24) are those given by Watson 
(loc. cit.), p. 441, 13.72, (2), and p. 150, 5.43, (1) respectively. 


5. The Nicholson-Watson integrals 
Thus far the analysis of the equations (7), (8) has been based on 
equation (10) satisfied by the algebraic function which I have called 
F,. I now go on to an analysis based on equation (9) satisfied by 
the analogous algebraic function £,,. 
Consider gee (7), and analyse it into the two equations 
{82—(u—v)*}h(z) = 42442) (z) 
oe 64 1)(5+-2)p(z). 
Taking ¢ = os y(2z), b= B..(e) ) we get the Mellin integral 


[elt 


The terms at the limits are 

40°{(1+-#)(8;+8,+-:3)—HK,,_,(2¢) £,,.,(t). 
These vanish as t> +0, if R(z) > 0, and as t> + 0, if |R(u—yv)| <1. 
We have therefore the solution of (7) 


E Ze) dt sic 
0 


if R(z) > 0, |R(u—v)| <1. 

If we write ¢ = zr in (27) and remember that every H(z) is of order 
z-1 at z= 0, we see that the integral (27) is itself of the order z-!. 
More conveniently we write t= cosechu and replace (27) by the 
integral ss 
| K,,-,(2z sinh uje“#*™ du. (28) 
0 

Now the general solution of (7) is a sum of constant multiples of 
the products 

J,4,; JY,» Yd, a 
At infinity the dominant terms of J, Y,, are 


i 
J,~,|(=)eose—ba—I9) Y »~,/(S) sions dow — 40. 
Tz :) 
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Thus only in ee J,I,+¥_Y, (29) 
and their linear multiples are the dominant terms free of the oscil- 
latory factors cosz, sinz. They have each the same order as (28), 
namely, 2-1, and therefore (28) is a sum of constant multiples of the 
two expressions (29). We can distinguish between the two expressions, 
and determine the appropriate constant multipliers by considering 
the second terms in the asymptotic expansions. We finally get that 


 y-y(22 sinh we“ du = }a*cosec(u—v)n(J,¥,—Jd,¥,,), (30) 


if R(z) >0, |R(u—yv)| <1. 

By changing the signs of pw, v and expressing J_,, Y_,, ete., in 
terms of J Yw etc., we deduce that 
( K,,-,(22 sinh u)cosh(u+-v)u du 
0 = (J.J, +Y,Y,+tan }(u—v)a(J,¥,—J,¥,)}, (31) 
subject to the same conditions on z, u—v. When p = v, these formulae 


become 


wo 


| K,(2zsinh wen ve du = - te ¥, 2 


ee 


0 


| K,(2zsinh u)cosh 2vu du = $7*(J?+-Y?), (33) 
0 
provided that R(z) >0 
Of these formulae (3¢ »), (31) are effectively those given by Watson 
(loc. cit.) as p. 445, 13.73 (4), (5), while (32), (33) are, of course, the 
formulae p. 444, 13.73 (1), (2). 


6. An integral involving W,.,, 

There are two other modes of splitting up equations (7), (8) which 
lead to results apparently new but not, I think, otherwise of interest. 
In the first place we analyse (8) into 

{8?—(u—v)?}p(z) = —42°(5+-2)4(2), (34) 
{(82—(u+-v)*Wle) = —2(8-+ Ihe). (35) 
These are equations of confluent-hypergeometric type. In pointof fact 
(34) has the solutions ete, yun ( 92) (36) 
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in Whittaker’s notation,* while solutions of (35) are 


eT sy(2*), e * Kyun?) re 


I am doubtful of the interest of these functions and I accordingly 
give the analysis only for the case » = v. We then replace equation 
(8) by the third-order equation (8B), which we analyse into 


5h(z) = —2°(2), 
(8°?—4v?)xb(z) = —427(8+ 1)h(z). 
Taking $(z) = e-!*, b(z) = e-**K,, (22), 


we consider the integral 


The terms at the limits are now 
£7(5,+8,+ 407+ 2)e2?"K (0), 


which, in virtue of the exponential factors, vanish both at = 0 and 
at t= 0, if R(z*)>0. We therefore have, as a solution of (8B), the 
integral a 


e-3et-* °K (#) dt 


t 


> 


0 
or, more conveniently, by writing /? = }zt’ and then dropping accents, 
the integral . 


[ e-aneren (get) [R(e) > 0), 


0 


(38) 


At R(z) = +0 the dominant term is 


J) [ « st-neg—¥ dt, 
Zi 0 


If we fold this integral over at t= 1 and write u = t!—t-!, it becomes 


io 4) 


2 /(=) [ « (u*+2)2 day, 


ny 


0 
which reduces without difficulty to 
4 —lp—2z 
wa*e®, 


Now the general solution of (8B) is a sum of constant multiples of 


* Whittaker and Watson, Modern Analysis, p. 336, 16.1. 
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I?, I,K,, K?. Clearly K? itself is the only combination which has 
the appropriate form at infinity, and comparing coefficients we write 


2K2(z) = [ e-WtK (zt) “J (39)* 


. 


if R(z) >0. 5 
The more general result, deducible from equation (8) when p +1, 
appears to be “ 
2K (z)K,(z) = [ eee, in -ehtt)K yo in. 
B v . »(u—v (utv)\2 v(2t)’ 
where R(z) > 0. , 
The corresponding analysis of (7) or (7B) does not seem to yield 
a rectilinear integral. 


7. An integral involving a hypergeometric function 
My last analysis is that of equation (7) into the equations 
(8—p+v)(5—p—v)h(z) = —4274(z), (41) 
(3+ p—v)(3-+-p-+v)p(z) = 2°(8+-1)(8-+ 22). (42) 
The first of these has the solution 
$(z) = 2J,(22); (43) 
the second has the solution 
ob(z) = 2-4) F(4— hu — dy, 1—du—4v; 1—7; 22). (44) 
Here again I give the analysis only for the simpler case in which 
y.=v. I then analyse (78) into the equations 
8(3—2v)h(z) = —42°(2), 
(5+ 2v)yb(z) = 2°(8-+1 (2), 
and take $(z) = 2” J,(2z), ab(z) = 2-2*(1—z?)"-4, 


so that we have the Mellin integral 


5 ace cle 
2) J,( = )eee—ay ©. 
[Jere sy 


Write ¢ = cosec wu, and consider instead the integral 


2” | J,(2zsin u)sin’u cos®’u du. 


The terms at the limits are now simply 
£41), (20)-"(1 274 


oc 2”+28in”+1y cos?” +1uJ, (22 sin wu). 


* A particular case of Watson, loc. cit., p. 439, 13.71, with Z =z, v = 2/t. 
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Near sinu = 0, this is comparable with sin?”*!w, and therefore the 
terms at the limits vanish when these limits are taken to be 0, 42, 
provided that R(v) > —4. We thus have as a solution of (7B) the 
integral ad 

2 f J,(2z sin u)sin’u cos”’u du, (45) 





0 

if R(v) > —}. 
If we expand J, under the sign of integration, we express (45) as 

a series of powers z?”+?", where r is a positive integer. The solutions 
of (7B) are J?, J, J_,, J?,, and of these J? alone, if v is not an integer, 
is composed of just the right powers of z. We therefore identify (45) 
with a numerical multiple of J? and compare coefficients. The leading 
coefficient in the integral is 

bar 

{T'(v+1)}+ | sin?’wcos*’u du = er ay? 
ar OTL (2y-+1)’ 
while that of J? is’ 2-2+4T'(v+ 1)}-. 
We therefore complete the identification in the form 
4a 
p 22” J,(22 sin u)si vu d (46) 
J *(z) = ——___— J,(2z sin w)sin’u COs*"u aU, ) 
( ) D(v+4)va ah 
0 
where R(v) > 3. 
If we similarly analyse (8B) into 

5(5— 2v)h(z) = 4274(z) 

(5+ 2v)yp(z) = 2°(5+ 1 )p(z), 
we merely write J for J everywhere in the preceding argument and 


so finally achieve the result 


har 
Dov 
[2(z) = ———__ T,(2z sin u)sin’u cos**u du, (47) 
D(v+3) 


0 
where R(v) > . 
These formulae, (46), (47), can be verified by expansion and direct 
integration. 
I have to thank Professor A. L. Dixon for references and sug- 


gestions. 











THE VALIDITY OF LAGRANGE’S EXPANSION 
By T. W. CHAUNDY 
[Received 1 October 1930] 


In the theory of the real variable Lagrange’s expansion is usually 
given in the approximative form that (under certain conditions) 


fy) fer+ > F gill KF (a)}+0 (a), (1) 


where y is defined implicitly by the relation 

y= a+xdg(y), (2) 
but as a rule nothing is said about the conditions, if any, in which 
the exact equation 


n 


fy) =fa+ >= A g@k@ (3) 
r=1 ~ 


is true. 

The proof of (1) consists in showing that the series on the right is 
the Maclaurin-Taylor expansion of f(y), regarded, in virtue of (2), as 
a function of x. To establish (3) we have therefore to prove that f(y) 
is exactly represented by this particular Taylor series. Conditions 
under which a function of a real variable is exactly the sum of one 
of its Taylor’s series have been published by Pringsheim* and by 
Young.t+ The form of condition which I shall use here is one which 


[ have given elsewhere. { 


The Taylor series eS 
satel n 
> Fay a) ; 
= n! 


if it converges, converges to F(z), aptcony that x,a both lie in an 
interval in which F(x) has derivatives F(a) of every order, and 

| FO (x)/n! [ln 
regarded as a function of x and positive integral n, is bounded. 


* Math. Annalen, vol. 44, p. 57. 

+ Quart. J. of Math., vol. 40, p. 157: ‘The fundamental theorems of the 
differential calculus’ (1910), p. 57. 

{ Messenger of Math., vol. 45 (1915), p. 115. 8S. Bernstein (Math.Annalen, 
vol. 75, p. 449) quotes this particular form as from Pringsheim ‘ Math. 
Annalen, vol. 45’ where I cannot trace it. Pringsheim (loc. cit.) gives a 
variety of equivalent forms, but, if I follow him, does not explicitly mention 
this particular form. 
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It is convenient to borrow from the theory of the complex variable 
and say that F(a) is ‘analytic’ in an interval in which every F(x) 
exists and | F'(a)/n!|/" is bounded. We are then reminded that the 
theory of the real analytic function is implicit in that of the complex 
analytic function and novelties are not to be looked for. There are 
reasons, however, for discussing the real analytic function indepen- 
dently of the more inclusive theory. 

For our present purpose, namely to prove the validity of the 
Maclaurin-Taylor expansion of F(x) = f(y), we have to prove F(z) 
analytic in the neighbourhood of x=0. We begin by proving in- 
ductively that the nth derivative F(x) is given by the formula 
= i é\"*fwlew)|" (4) 

(1—a"(y) 6y) = 1—ad'(y) ” 
presupposing the necessary derivatives on the right. 


F(x) 


Now from (2) we have 
da Y d 0 0 0 
= a a me 2a ae oy 
Writing P|Q for the alternant PQ—QP of operators P, Q, we have 
14 2 6 8 
éx|1—ad’ dy (1—ad’)? ey’ 


$ | 1 eo -—# @ 
1—ad’ éy|1—ad’ ey (1—ad")? dy’ 
d ] a 
dx|l—ad’ dy” 
i.e. these two operators are commutative. Thus for operation on any 
differentiable function %(y) we get 
dj 1 e\" wy) | 1 d\*td 
dx\1—ad’ ey) 1—axd’ |1—ad’ ey} eh) 
_f{ 1 @\""* oy) 
\l—ad’ dy! 1—ad’’ 
Hence, if we here write ¢’(y) = f’(y)[¢(y)|", which is legitimate, since 
i'¢", being differentiable, is also integrable, we have proved that 
d 
—Y,, = Vn 
dx n n+l 
where Y,, stands for the expression on the right of (4). Again 
, (y)fi(y) 4d 
Y,= $ wf F = f(y), 
l1—ad'(y) = da 


and so, by addition, 
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which is correct. The induction is therefore complete and the formula 
(4) proved. 

We use (4) to prove the essential theorem: 
(1) If for some x,y 
fmy)|, Ay : <n!A B" (every n >0), 
and \l—ad'(y)| 1< C, 
then pe 1)| <n!A, Bt (every n > 0), 
where A, B, C, Ay, B, denote positive constants independent of x, y, n 
The proof of (1) is based on the two lemmas: 
) If |fo(y)|, Pry) <n!AB" (every n > 0), 
and p is integral, 
then |D"(f'h?)| <n! AP+1 Bi+1 (every n, p > 0), 
where D denotes differentiation in y, and A,, B, are positive constants 
independent of x, n, p; 
and 
3) If fm(y)| <n!A,B", \(y)| <n!AB" (everyn>0), 
\dz|- 


and |—| nag GC, 
dy 


anf) 
then CI <mn!A,A, Bt (every n > 0), 
where A,, B, are independent of x, n, Ag. 
To prove (2) we have by Leibniz’s formula 


pn cf — aii oe. —8) 


s=0 


D"(f'd)| <nlA2Brt > (s+1)= = A2 Bn, 


s=0 


D% f $2) = i> & (f’ ¢) aaa : 


s=0 


2 ee 
D(f'$2)| <n!A® Bu s = = " 3)! As Bri, 


s=0 
Proceeding inductively in this way we get at length for any positive 
integer p 
: » . (n+p+1)! 
\D'(f p?)| < = AP+ 1 Br 1 
| (p+1)! 
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+1 ; : : 
Now wee a is a coefficient in the expansion of (1+-2)"+?+1 and 
m!(p+ 


is therefore certainly less than 2”+”+1, the sum of the coefficients in 
the expansion. We thus have that 


|D"(f’p?)| <n!(2A)?+1(2B)” 
which proves the lemma (2). 
To prove (3) we write 
d” (dz —(n+s8) qd" . 
a > Lanza ae) dy a (9) 
r=0 s=0 y Y 
d*z dz : . 
is a function of arguments —~, —,... only. Differentia- 
dy” dy*® 
tion and comparison of coefficients give the recurrence-formula 
7, : | 2 
Mua 17,8 — Lara t Diufing-1.5 + (n i ci ] )(D3z)Zn»—1,6-1 (6) 
which, used inductively, shows Z,, ,, to be a polynomial with positive 
numerical coefficients. 
Let us write Um, = (Dirz)/m!, (7) 
so that Z,,,, becomes a polynomial in u,, (m = 2, 3,...) with positive 
coefficients. Its degree and weight in w,, being its orders in z and 


where Z 


n,?,8 


y—1 respectively, we see from (5) that Z,,,., is homogeneous of degree 
s and isobaric of weight r+-s. 


Since, from (7), Dyin = (M+ 1 Jina 


the sum of the coefficients in D, [lu is 
> (r+ la, =i+u, 
where i, w are the degree and weight of Ilu%. It follows that, if 
c is the sum of the coefficients in some f(g, u3,...) which is homo- 
geneous and isobaric of degree 7 and weight w, then (i+w)ec is the 
sum of the coefficients in D,f. 
Hence, if c,,,, is the sum of the coefficients in Z,,,.,, the sum in 


ns 
D,Zn7,3 18 (T+-28)¢,,y,5- From (6) we therefore derive the recurrence- 
formula 


Cn+t18 — Crns »+(r+ 28—1)c ‘nf ast 2(n+s8—I1)ec 2,f—1,8—19 
the coefficients being additive, since they are all positive. This 


recurrence-formula, together with the boundary values 


Cnyr,s — 0, if r>n—1 or s>y7, Case = 1, 
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defines inductively a unique c,,,,. By scientific, rather than mathe- 
matical, induction I discover the formula 


i 
mrs (n—r—1)\(r—a)!a!(s—1)! 
which can be verified by sheer substitution in the defining equations. 
Since |u,,| <<AB” by the hypothesis of (3), we have 


mi 
_ = twee, 





while from gi we have 


dnf|_ ~(2n—1) — | qn- % \dz|"-r-1 J 
- Z By poset 
an > dy” + dy > ns y 


s=0 
Since also by i ~~ viothesis of (3) 
-1 
a"S| < m\A,B™, <o, |@|<as, 
idy™ dy dy 
we get 
df w n—1 r ; : a 
; 1 |< (2n-1 » 2 (n—r)!A,B"-(AB)"7-1 > Cu rg4*°Brt4(A By 8 
az" r=0 s=0 -° 
-1 
ie.  <A,A"-(BC)"-'5 (n— 
r=0 


Now 


n—1 





= (n+s—1)\(r—1)! 
> (n—n)! >e nrs = > (n—1)\(n—r ) >. ja—Tyalr—ay \(s—1)! 


r=0 


This is the pee me is t”-1 in 
n—-1 


2 (n—1)!(n—r) > aie ay (1-+4)"+8-1 
s=0 : : 


r=0 
n—-1 


=> (n—Il(n—r)(1+4(2+ty 
r=0 
This expression is less, term by term, than 
nl(L-+t)e'S (2+ty 3 = al(1--4)-4((2-+4)"-4— 1. 

The coefficient of ¢”-1 is certainly less than the sum of the coefficients, 
and so less than n!2”-13"-1, Hence finally 

os \< n!Ag! (6.4 B2C?)” 

2"  6ABC ’ 
which proves the ie (3). 
To prove (1), write as before 


v'(y) =f Wiley? 
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for positive integral p. Then cad (2) 
Wry) | < (n—1)!AP+1 BY. 
Write 2(y) = Sa as 
where x is now to be rh as a constant. Then 
-r"(y) | - |P™(y)| <n!A, By, 

and therefore by (3) 

(a 

| 1—2d'(y) ey 
Now take p = n and we have 


(n) (> _|f 1 c I; 
I (=) =|) 7 xd) dy! bly on 


< m!AP+1A, BY. 


which proves (1). 
It remains to translate (1) into its permanent form: 
If the functions f(y), d(y) are analytic near y=a, and if y(x) 
is defined implicitly by the relation 
y=—a+-xd(y), 
where, if y(a) is many-valued, we mean that branch which has the value 
a at x= 0, then y, or more generally f(y), is an analytic function of x 
near x = and is exactly equal to the sum-to-infinity of its Lagrange 
series. 
Since ¢’(y) exists and is continuous near y = a,and 1—2x¢'(y) does not 
0, then, by the theory of implicit functions,* the relation 
= a-+-ad(y) 
defines near 2 = 0 a one-valued, continuous function y(#) such that 
y(a) = 0. If x is sufficiently small, we can therefore secure both that 


1—a¢'(y) > some 1—C>0 


vanish at x 


and that y(a) lies in a prescribed neighbourhood of y= a. The con- 
ditions of (1) are then satisfied and accordingly F(x) is analytic in 
this neighbourhood of «= 0 and is represented accurately by its 
Maclaurin-Taylor expansion, which is Lagrange’s series for f(y) 

The foregoing analysis can be extended without difficulty to show 


that certain mixed functions of x,y such, for instance, as 

[1—ad'(y) |" 
are also analytic functions of x near x 
the appropriate modification of Lagrange’s series. 


-Q and are represented by 


* Cf. Goursat, Cours d’analyse (1910), vol. i, § 34, p. 81. 











